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I.  Introduction 

In  this  review,  we  will  examine  the  possibility  of  utilizing  energecic  ton  beams  as 
cools  to  characterize  the  chemistry  and  structure  of  solids  and  surfaces.  This  approach  may 
appear  somewhat  unusual  to  most  spectroecoplsts  who  art  generally  concerned  with  the  iuantuo 
mechanical  Interaction  of  electromagnetic  radiation  with  aatcer.  But,  chere  have  been  a 
number  of  developments  that  hava  surfaced  over  the  late  few  years  which  have  attracted  our 
attention  to  this  particular  field  and  which  merit  close  consideration. 


The  Idea  behind  Che  experiment  of  interest  is  illustrated  in  Figure  l .  We  begin  with  a 


Figure  l  -  Positions  of  the  atoms,  (a)  Before  the  primary  ion  (the  lone 

sphere  above  the  solid)  has  struck,  (b)  Consequences  of  a 

single  ion  impact.  The  two  atoms  ejecting  to  the  left  form 
a  dimer.  For  graphical  clarity,  only  a  selected  group  of  atoms 
is  shown,  and  their  size  la  arbitrarily  drawn  for  best  graphical 
clarity. 

single  crystal  substrate  represented  by  an  ordered  array  of  spheres  containing,  in  this 
case,  chree  layers.  The  atom  above  the  crystal  in  Figure  la  represents  the  Incident 
particle  whose  kinetic  energy  is  generally  bacumen  a  few  hundred  and  a  few  thousand  electron 
volts.  After  the  collision,  the  atoms  In  the  crystal  begin  to  move  in  various  directions  as 

the  energy  of  the  incident  particle  is  shared  between  these  atoms.  In  Figure  lb,  after  less 

than  10-12  s,  some  particles  have  ejected  into  the  vacuum  and  may  therefore,  in  principle, 
be  detected  in  the  laboratory. 

The  question  is  -  can  vm  learn  something  about  the  solid  by  a  detailed  analysis  of  the 
motion  subsequent  to  the  primary  bombardment  event1?  The  answer  and  implications  of  such  a 
question  are  certainly  not  obvious.  In  fact,  to  illustrate  the  difficulty  it  is  appropriate 
to  make  an  analogy  to  the  game  of  pool.  Bare,  a  cue  ball  is  Incident  upon  a  triangular 
array  of  IS  target  spheres  placed  on  an  isolated  surface.  What  we  would  like  to  know  is  - 
suppose  a  player  comes  into  the  room  with  Che  pool  cable  and  then  initiates  the  action. 

After  all  the  balls  scop  moving,  he  leaves  the  room  and  a  second  observer  enters  who  is 
asked  to  reconstruct  the  original  configuration  of  Che  balls.  He  is  asked  not  only  to  reveal 
their  original  structure  but  also  to  place  the  numbered  balls  in  their  respective  sequence. 
Thue,  was  Che  eight  ball  originally  touching  tha  two  ball?  If  we  could  do  this  on  an  atomic 
scale,  we  could  obviously  learn  a  great  deal  about  the  structure  and  interaction  forces 
between  atoms  near  the  point  of  impact  of  the  bombarding  particle. 

This  general  subject  has  been  Important  outside  of  the  field  of  chemistry  for  nearly  a 
century.  The  epplications  appear  in  such  diverse  fields  as  ion  implantation,  particularly 
for  doping  seml-conductora ,  fuaion  reactor  design  and  the  determination  of  the  effect  of 
well  reactions  on  plasma  temperatures,  and  astrophysics  and  the  study  of  the  influence  of 
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Che  solar  wind  (He'*')  on  extraterrestrial  surfaces.  The  lnceresc  may  be  focussed  on  the 
damage  Imparted  to  Che  solid  by  the  Incident  particle  or  to  the  analysis  of  the  ejected 
material. 

There  Is  a  long  history  Co  the  development  of  Che  current  state  of  this  field.  Before 
1970  most  of  che  research  was  carried  out  by  physicists  who  really  mere  not  concerned  with 
any  of  che  possible  chemical  Implications.  Furthermore,  aosc  of  the  experiments  were 
completed  under  rather  Ill-defined  conditions.  Thus,  we  have  chosen  not  to  comprehensively 
review  all  published  literature  but  to  cry  to  congeal  from  many  fragmentary  pieces  of  data  a 
scenerlo  which  should  be  of  interest  in  che  development  of  new  approaches  to  the  chemical 
and  structural  analysis  of  a  wide  variety  of  materials. 

There  are  a  number  of  excellent  reviews  of  related  copies  for  those  readers  Interested 
in  filling  In  Che  gaps.  First,  there  Is  a  landmark  publication  by  Carter  and  Colligan  (1) 
which  Is  comprehensive  and  complete  through  1967 .  !iext,  there  have  been  several 
publications  Chat  have  appeared  during  the  last  5  years  which  are  reprints  from  major 
symposia.  These  Include  an  International  Workshop  on  Inelastic  Ion-Surface  Collisions  held 
in  July  of  1976  at  Bell  Laboratories  In  Murray  Hill,  Mew  Jersey  (2),  Che  proceedings  of  the 
8th  International  Conference  on  Atomic  Collisions  In  Solids  held  In  August  of  1979  in 
Hamilton,  Canada  (3),  che  Second  International  Conference  on  SIMS  In  August  of  1979  at 
Stanford  University  (4),  an  International  Symposium  on  Sputtering  held  in  Vienna  In  April  of 
1980  (5)  and  an  International  Workshop  on  Inelastic  Ion-Surface  Collisions  held  in 
Feldkirchen-Westerham,  Germany  in  September  of  1980  (6).  Finally,  there  are  several  journals 
usually  not  familiar  co  Chose  outside  of  the  field  Including  Radiation  Effects  and  che 
Journal  of  Nuclear  Instruments  and  Methods  where  a  number  of  relevant  articles  may  be 
found. 

A  further  point  Is  chat  we  shall  ignore  completely  che  ultlmace  fate  of  the  lncidenc 
particle.  This  aspect  of  the  problem,  however,  also  provides  complementary  Information  and 
there  are  fascinating  experiments  proceeding  along  these  lines.  For  example,  che  techniques 
of  Ion  scattering  spectrometry  (6),  MsV  backscattering  (7),  and  atomic  beam  diffraction  (8) 
are  providing  intimate  details  about  che  structure  of  solid  surfaces. 

The  development  of  this  article  will  proceed  along  the  following  path.  We  will  begin  by 
citing  a  few  of  che  key  experiments  which  Influenced  most  strongly  che  currenc  direction  of 
che  field.  Nexc,  we  shall  review  some  of  the  major  theoretical  developments  that  have  led 
to  che  current  detailed  understanding  of  what  happens  during  the  bombardment  event,  and 
finally  wt  will  elucidate  soma  of  che  applications  chat  should  be  of  lnceresc  to  chemists  Lr 
ehelr  conetanc  search  for  new  types  of  analysis  mechods. 


II.  Ion  Bombardment  Methods 


To  characterize  Che  ejection  or  material  from  the  solid.  It  la  possible  to  devise 
experimental  techniques  to  detect  both  the  neutral  particles  and  those  thee  are  Ionized 
directly  by  the  surface.  Although  It  la  easier  to  detect  snail  quantities  of  ejected  Iona, 
the  neucrals  have  more  Interest  froa  a  theoretical  point  of  view.  In  this  section,  the 
experimental  approaches  to  characterizing  the  ejected  neutrals  and  ions  will  be  presented. 
The  emphasis  will  be  placed  on  the  use  of  single  crystal  samples  since  che  quesclon  of 
atonic  structure  la  central  to  this  article. 

A.  Sputtering 

As  shown  in  figure  1  when  a  keV  ion  strikes  a  solid,  there  is  sufficient  momentum 
transferred  to  the  target  atoms  that  considerable  nuclear  motion  la  Initiated.  A  fraction 
of  che  moving  atoms  may  obtain  a  sufficient  component  of  momentum  oriented  Into  the  vacuum 
that  they  are  able  to  overcome  che  surface  binding  force  and  ejecc  from  che  solid.  Some  of 
these  particles  emerge  directly  as  positive  or  negative  ions,  while  others  are  ejected  as 
either  ground  or  excited  scace  neutral  atom a  as  illustrated  In  Figure  1.  The  emission 


Figure  2  -  Schematic  representation  of  the  Ion  bombardment  experiment. 

The  primary  particle,  ?,  generally  Is  given  a  few  hundred  to  a 
few  thousand  eV  of  kinetic  energy. 

process  and  the  associated  radiation  damage  In  che  solid  are  often  referred  to  as 
spuccerlng,  a  term  left  over  froa  che  earliest  observations  of  the  phenomenon  bsck  In  che 
1850's.  These  workers  were  concerned  about  the  erosion  of  cathodes  in  gas  discharge  sources 
(9).  According  to  a  recent  review  CIO),  the  word  "sputtering-  intimates  a  microscopic 
appreciation  for  che  ejection  process  sines  it  implies  chat  matter  it  emitted  in  discrete 
blobs  from  the  areas  of  the  surface  that  are  struck  by  the  bombarding  Ion.  In  the  modern 
sens*1,  the  dictionary  definition  of  the  word  is  "  To  emit  sallvs  from  the  mouth  In  small 
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particles  hardly  a  scientific  Cera.  It  appears,  however,  chac  the  old  terminology  Is 

desdned  co  remain  with  us  for  soma  doe.  Historically,  then,  the  first  studies  of  Ion 
bombarded  solids  Involved  the  investigation  of  the  sputtered  species. 

Moat  Ion  bombardment  scudles  up  until  the  mid  1960's  are  spent  collecting  seemingly 
unrelated  determinations  of  the  sputtering  yield  as  a  function  of  all  possible  variables. 
Here  the  yield,  Y,  la  defined  as 


Y-S/i0  (1) 

where  S  Is  the  number  of  ejected  particles  (atoms  plus  ions)  from  the  target  and  10  Is  the 
nunber  of  Incident  Ions .  Y  was  found  co  vary  with  the  angle  of  Incidence  of  the  primary 
ion,  the  atomic  weights  of  all  Involved  species,  the  primary  Ion  energy,  the  primary  Ion 
current  and  the  angle  of  che  ejecting  particles  (1).  The  value  of  Y  was  nearly  Independent 
of  temperature  at  temperature  values  wall-below  the  melting  point  of  che  target  (1). 

The  investigations  mainly  Involved  measurement  of  tha  weight  loss  of  the  target.  Other 
analysis  methods  Included  measurement  of  changes  In  the  resonant  frequency  of  a  crystal  on 
which  the  target  Is  plated  (11),  noting  changes  of  the  target  resistance  (12)  and  use  of 
radiotracers  (13).  Reproducible  results  were  extremely  difficult  to  obtain  since  vacuums 
were  only  ~10-^  corr  and  impurities  could  have  a  large  influence  on  the  results .  In 
addition,  the  morphology  of  che  sample  Influenced  Y  and  this  property  would  generally  change 
during  che  measurement.  Nevertheless,  there  were  several  key  experiments  that  have 
Influenced  present  day  concepts  of  ion/ solid  interactions. 

The  Influence  of  Y  on  the  surface  morphology  provided  che  first  clue  to  the  fact  that 
the  momentum  dissipation  of  the  primary  Ion  was  sensitive  to  the  atomic  structure.  In 
Figure  3,  for  example,  Is  shown  an  etch  pattern  of  s  Cu(ll,3,l)  target  and  a  single  crystal 


Fljura_3  -  Electron  micrograph  of  surfaces  bombarded  by  9  keV  Ar+. 

(a)  Cu( 11,3,1),  magnification  2000x  (b)  SI  crystal, 

magnification  SOOOx.  From  reference  14. 


silicon  surface  bombarded  by  9  keV  Ar+  ions  (14).  In  both  cases,  the  notable  feature  is  one 
of  facetting  which,  for  Cu,  seems  to  emerge  from  the  base  of  etch  pits.  The  theory  that 
different  facets  of  a  polycrystalline  surface  exhibited  different  yields  was  proposed  by 
Uehner  some  years  earlier  (IS, 16).  In  fact  Uehner  found  preferential  etching  patterns  which 
reflected  the  atomic  symmetry  of  the  target  crystal  face.  Since  these  studies  pictures  of 
many  spectacular  surface  features  have  been  published,  including  those  of  great  columns  and 
cones,  like  chose  of  Figure  1  reminiscent  of  the  surfaces  of  distant  worlds.  These  changes, 
however,  indicated  that  the  observable  could  be  a  function  of  primary  Ion  dose  (dose  • 
ions-sec'^-cm~l  x  else).  Ic  became  apparent,  then,  chat  techniques  which  required  the 
smallest  dose  would  be  most  desirable. 

Ulth  Improvements  In  vacuums  and  primary  Ion  sources,  th*  quality  of  yield  measurements 
began  co  Improve  markedly.  Initially,  it  was  found  Chat  T  generally  increased  with  the 
kinetic  energy  of  che  primary  loo  to  a  maximum  value.  At  even  higher  klneclc  energies,  the 
yield  then  began  co  decrease,  presumably  since  che  primary  ion  penetrated  too  far  Into  the 
target  to  cause  ejection  of  surface  atoms.  Of  particular  Interest  was  that  there  is  a 
scrong  dependence  of  yield  on  che  orientation  of  single  crystal  samples.  In  Figure  4,  for 
example,  is  presented  the  yield  vs.  primary  kinetic  energy  curve  for  Ar*  normally  Incident 
on  CuUll),  Cud  10),  Cu(100)  and  polycrystalline  Cu  (17-19).  yote  that  at  least  above  l  keV 
Ar*  Ion  energy,  T(lll)  >  T(IOO)  >  T(llO)  and  chat  the  maximum  yield  can  approach  10  ejected 
Cu  atoms  per  incident  Ar+  Ion.  This  same  trend  has  been  reported  for  Ag  (19)  and  Au  (20) 
indicating  the  Importance  of  crystallographic  effects.  Similar  variations  of  yield  with 
angle  of  Incidence  have  been  reported  as  shown  In  Figure  5  (21)  where  the  yield  minima 
correspond  to  open  crystallographic  directions.  We  must  admit  that  these  results  are 
somewhat  surprising,  especially  In  view  of  the  rather  scrong  effects  che  lncidenc  Ion  may 
have  on  che  surface  morphology  as  demonstrated  in  Figures  1  and  3.  It  would  be  particularly 
valuable  to  repeet  these  studies  where  surface  dsmsge  effects  do  not  occur. 

Crystallography,  however,  clearly  is  an  important  factor  in  sputtering  which  may  be  the 
reaaoa  chat  moat  workers  during  che  1960's  spent  much  effort  analyzing  single  crystals. 
Probably  che  most  spectacular  discovery  of  chat  ere  was  chat  of  Wehner.  Using  the  very 
simple  discharge  tube  apparatus  shown  in  Figure  6  (22),  he  utilized  a  Hg+  ion  discharge  to 
bombard  a  Cu  crystal  and  collect  che  sputtered  atoms  on  a  flat  plate.  An  example  of  the 
type  of  pattern  that  could  be  collected  using  this  procedure  is  given  in  Figure  7  (IS, 23). 
The  patterns  showed  clear  evidence  for  large  angular  anisotropies  In  the  yields  which  also 
reflected  the  synsstry  of  the  atoms  in  the  surface  layer.  Since  these  papers,  more  then  100 
studies  have  appeared  in  the  literature,  detailing  che  effects  for  e  wide  variety  of 
crystals  under  many  ion  bombardment  conditions.  These  studies  have  waned  in  recent  years 
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Figure  7  -  Angular  distributions  of  atoms  ejected  from  CulOOl)  from 
reference  17. 

due  to  che  lack  of  a  cohesive  theory  to  explain  the  anisotropy  and  due  to  the  experimental 
difficulties  of  maintaining  an  unconcamlnated  and  undamaged  crystal  surface. 

An  Improved  approach  has  recently  been  developed  for  studying  these  anisotropies. 
Specifically,  the  sputtered  material  Is  collected  on  an  aluminum  place  and  subsequently 
analyzed  using  Rutherford  backscacterlng ,  a  highly  quantitative  and  sensitive  measurement 
technique  which  requires  nearly  two  orders  of  magnitude  smaller  primary  ion  doses  than 
previous  methods  (24).  With  chls  apparatus.  It  has  been  found  that  for  gold  (111) 
bombarded  by  4  keV  s»*  chac  more  chan  302  of  Che  spuccered  material  Is  emltced  Into  Che 
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spots.  The  shapes  of  the  spots  ere  nearly  independent  of  the  projectile  use  and  energy  for 
energies  greater  than  ~6  keV  (26).  These  findings  contrast  previous  assertions  that  these 
anisotropies  constitute  only  a  minor  part  of  the  sputtering  process  and  may  therefore  be 
Ignored  (10). 

An  additional,  important  place  of  information  which  provides  insight  into  Che 
sputtering  process  is  the  velocity  distribution  of  the  ejected  atoss.  These  were  first 
reported  by  Thompson  and  Nelson  (13)  using  a  time  of  flight  analysis  of  Che  atom  velocities 
and  later  by  Stuart  and  Wehner  using  a  fluorescence  net hod  ( 16) .  The  distributions  as  shown 
in  Figure  8  are  similar  for  most  uterlals  in  that  the  curve  is  0  at  0.0  eV,  peaks  s'.  Che 
fairly  low  energy  of  2-5  eV  and  then  tails  off  with  an  energy  dependence  of  somewhere 
between  1/E  and  l/E^.  The  total  yield,  T,  then  can  be  expressed  in  terms  of  an  energy 
spectrum  where 

T»t*y<E)dE  (2) 

0 

where  y(E)dE  is  Che  differencial  secondary  particle  yield,  and  E^  is  Che  oaxitnuo  energy  th&c 
any  of  the  ejecced  particles  has  attained.  The  energy  distributions  are  markedly  different 
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Figure  8  -  Energy  distribution  of  atoms  ejected  in  the  [1101  direction 

from  Che  (110)  surface  under  bombardment  by  norully  incident 
Kr  ions.  From  reference  16. 
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from  chose  expected  for  thermal  evaporation  processes  where  Che  full  wldch  at  half  maximum 
of  che  peak  would  be  Maxwell-Boltzman  in  shape  and  considerably  smaller. 

Although  these  studies  have  provided  much  of  the  fundamental  Information  that  we  have 
regarding  che  ion  bombardment  event,  the  scope  of  Che  techniques  la  severely  limited  by 
their  ultimate  sensitivity  and  general  lack  of  chemical  specificity.  A  number  of  approaches 
have  been  tested  In  order  to  obtain  better  Information  from  the  ejected  neutrals.  One  Idea 
Involves  che  post-lonizacion  of  che  ejected  species  either  by  eleccroa  bombardment  (25)  or 
bv  a  low  energy  plasma  discharge  (26).  The  Ion  which  Is  formed  Is  Chen  detectable  vlch 
single  particle  detection  limits  using  a  mass  spectrometer.  This  concept  has  been  exploited 
by  Konncn  et  al.  (27)  who  employed  a  time  of  flight  technique  to  determine  particle 
velocities,  but  tries  utilized  an  ionization  scheme  to  detect  the  mass-selected  beam.  This 
procedure  worked  well  for  alkali  metals  but  did  not  provide  sufficient  sensitivity  for  ocher 
particles  (28,29).  In  fact  most  poat-lonizacion  methods  are  hard  pressed  to  compete  with 
the  older  sputtering  methods  based  on  weight  loss  analysis.  The  efficiency  of  Ionization  by 
electrons,  for  example.  Is  generally  only  about  10~*.  In  addition  the  Ionization 
cross-section  is  strongly  dependent  on  the  particle  velocity,  a  severe  difficulty  for  these 
experiments  because  of  the  broad  energy  distribution  of  the  ejected  particles.  And  finally, 
the  process  of  electron  bombardment  may  in  Itself  cause  fragmentation  of  ejected  molecular 
clusters,  muddling  the  picture  still  further. 

The  use  of  a  hoc  maxwelllan  electron  gas  maintained  by  electron  cyclotron  wave 
resonance  has  shown  considerable  promise  as  an  approach  to  study  the  sputtering  phenomenon 
(26),  Ionization  of  neutrals  occurs  with  an  efficiency  approaching  IX.  In  addition, 
sputtered  neutral  mass  spectrometry  (SUMS)  has  been  successful  In  measuring  dimer  to  monomer 
ratios  for  a  number  of  different  metals.  For  example,  Oechsner  found  that  the  Cuj/Cu  raclo 
from  polycrystalline  Cu  was  approximately  10Z  (30). 

Perhaps  one  of  the  most  promising  methods  for  analysis  of  the  neutral  particles  Is  via 
che  use  of  laser  spectroscopy.  One  such  approach  Involves  the  excitation  of  the  ejected 
atom  by  a  laser,  followed  by  che  measurement  of  che  doppler-shif ted  fluorescent  Intensity 
(31).  Using  this  approach,  very  high  quality  energy  distributions  can  be  obtained  for 
sputtered  neutrals  using  Incident  Ion  doses  small  enough  Co  svold  significantly  altering  che 
surface.  The  block  diagram  for  a  Cyplcal  apparatus  Is  shown  In  Figure  9,  with  che  resulting 
energy  distribution  shown  In  Figure  10  (32).  It  Is  likely  that  with  che  advent  of 
ultra-high  vacuum  Instrumentation  and  the  evolution  of  sensitive  methods  of  analysis  of  the 
neutrals  that  very  high  quality  sputtering  dace  from  well-defined  solid  surfaces  will  soon 
be  available.  At  the  present  time,  however,  it  Is  not  possible  to  perform  yield 
measurements  wlthouc  damaging  the  crystal  surface. 
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Figure  9  -  Diagram  of  apparacua  used  in  the  doppler-shifted  laser 

spectral  experiments  used  to  measure  the  energy  distribution 
of  the  neutrals.  From  reference  32. 


Figure  10  -  Number  density  velocity  distribution  of  neutral  ground 
state  zirconium  atoms  sputtered  by  normally  Incident 
1.3  keV  .  The  (0)  represent  daconvoluted  date 
averaged  ovar  several  runs.  The  solid  curves  are 
obtained  from  Thompson's  squatlons  (13).  An  average  error 
Is  displayed  In  the  right-hand  corner.  From  reference 
32. 


3.  Secondary  Ion  Emission 


It  has  been  recognized  Cor  more  Chan  SO  years  (33)  that  a  fraction  of  Che  particles 
ejected  from  an  Ion  bombarded  target  are  either  positive  or  negative  Iona.  The  secondary 
positive  (negative)  ion  yield  rK-)  is  defined  as 

Y*  -  T  (3) 

where  ST*"  is  the  probability  chat  the  ejecting  neutral  species  will  leave  the  solid  In  sn 
Ionized  state.  A  similar  expression  could  be  written  for  negative  Ions.  It  Is  possible 
Chat  K+  is  a  function  of  energy  so  that  combining  Eq.  (2)  and  (3)  we  obtain 

Y*  -  fV  (E)  y(E)dE  .  (4) 

o 

For  many  systems,  the  value  of  y  may  be  large  enough  to  permit  detection  of  very  low 
concentrations  of  surface  spades.  Much  of  the  early  interest  In  secondary  Ion  emission  In 
fact  stemmed  from  this  potential  analytical  application.  Initially,  a  lot  of  effort  was 
spent  developing  che  ion  microscope  (34)  and  the  scanning  loo  microprobe  (35).  Here,  the 
secondary  Ions  were  directly  Imaged  into  a  double  focussing  mass  spectrometer  and  detected 
with  single  Ion  limits.  Evaluation  of  che  secondary  Ion  Intensity  T*  [counts/sec]  recorded 
during  bombardment  of  the  sample  Is  more  complex  than  for  the  neutral  species  as  determined 
by  weight  loss,  but  It  Is  easy  to  see  the  potential  sensitivity  of  this  approach.  For  this 
case  ignoring  angular  effects  (36) 

s2 

1+  •  l0rM8(^)T(E0,R )/  y+(E')dE'  (5) 

*1 

where  i0  is  che  flux  of  primary  ions,  yjj  is  the  isotope  abundance  of  element  M,  S(M*)  is  che 
detector  efficiency  which  depends  on  the  mass  and  velocity  of  the  ion  (37)  and  T(E0,R)  Is 
che  transmission  of  the  ion  optical  system  which  depends  upon  energy  poslton  E0+AE/2  and 
width  1E"E»-Sj  of  the  energy  window  and  che  mass  resolution  of  che  mass  spectrometer.  For 
typical  values  of  10  ■  10“^  amps  of  Ar*  Ions  spread  over  10”2  cm2  ( 6x10 1 3  ions-see_1— cm-2) , 
rjj2T  *  10"2,  E3  *  10  eV,  4E  »  5  eV  and  the  Integral  In  Eq.  5  equal  to  approximately  10"2, 
che  t*  «  6x10®  counts/sec.  With  a  detection  Halt  of  -1  cps.  It  is  clear  that  measurements 
In  che  parts  per  million  (ppm)  atomic  percent  range  can  be  completed.  In  addition,  there  Is 
also  sufficient  sensitivity  to  detect  Ions  ejected  from  a  very  small  area.  Early 
development  of  these  secondary  Ion  mass  spectrometric  (l.e.  SIMS)  methods,  however,  was 
hindered  mainly  due  Co  the  che  very  high  cost  of  che  Instrumentation. 
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Th«r«  are  two  additional  disadvantage!  of  this  approach.  First,  it  turns  ouc  to  ha 
vary  difficult  to  obtain  quantitative  ion  yield  measurements.  This  is  due  to  the  difficulty 
in  evaluating  ST  and  to  the  large  variations  found  in  y+(E)  in  Eq.  5  (or  K+  in  Eq.  3)  with 
the  saaple  matrix.  For  example,  Che  yield  of  Ml*  from  Ni  is  103  times  smaller  than  the 
yield  of  Hi*  from  S10  (38)  under  otherwise  identical  conditions.  To  understand  these 
results  it  is  necessary  to  have  a  theoretical  understanding  of  the  secondary  ion  formation 
mechanism.  As  we  shall  see  in  Section  III,  there  are  a  number  of  possible  approaches  to 
explaining  che  ionization  process,  although  none  have  yet  proven  to  be  very  satisfactory. 

A  second  major  problem  with  SIMS  as  presented  above  is  that  although  the  primary  ion 
current  is  small,  che  primary  ion  flux  is  extremely  large,  typically  10  to  100  uaap/em2,  or 
lO1"*  co  1016  incident  ions-sec-3.  If  we  assume  there  are  roughly  1033  acoma-cm”2  on  a 
typical  mecal  surface  and  che  sputtering  yield,  Y,  is  -1,  Chen  the  saaple  is  being  eroded 
away -at  the  rate  of  0.1  to  10  monolayers-sec~3 .  The  magnitude  of  che  sample  damage  Is 
clearly  quite  large  and  the  pr<^c><-  .  ;'  -r  performing  a  chemical  analysis  on  a  molecular 
syscem  would  appear  dim. 

During  che  lace  196o'<  ,  t  »..  ••'lea  was  ingeniously  overcome  by  Bennlnghoven  and  his 
co-workers  (39,10)  who  simp;.-  , i--’ Vised  to  expand  the  size  of  the  primary  ion  beam  and  reduce 
the  total  dose  of  lone  leas  t.  in  -10 13  atoms-ew"2.  By  using  incident  currents  of  ~ICr- 
to  10-30  amps-cm'2 ,  then,  he  could  study  the  chemical  nature  of  a  surface  for  minutes 
without  altering  che  surface  composition.  The  erosion  race  was  effectively  reduced  to 
nr5-ur4  monolayer-see"3.  In  addition,  since  a  wide  area  was  being  Imaged  into  che  mass 
spectrometer,  he  was  able  to  considerably  enhance  sensitivity  by  utilizing  a  quadrupole  mass 
filter  lnscead  of  a  magnetic  sector  (41).  This  implementation  was  also  Important  since 
the  energy  distribution  of  che  ejected  species  extends  over  quite  a  range  as  shown  in 
Figures  8  and  10  and  quadrupoles  are  generally  insensitive  to  ion  energy.  The  cost,  not 
incidentally,  la  considerably  lower  and  che  equipment  can  be  made  compatible  with  ultra-high 
vacuum  (UHV  or  10"30  torr)  conditions.  A  typical  experimental  set-up  (42)  is  shown  in  Figure 
11. 

Thus,  there  ha.'*  evolved  several  approaches  co  characterizing  che  phenomenon  of 
secondary  ion  emission.  There  is  che  small  beam,  ion  microprobe  SIMS  method  chac  provides 
mlcroanalysls  and  high  sensitivity  but  rapidly  erodes  che  sample  surface.  At  che  che  other 
extreme,  there  is  the  large  beam,  low  dose  SIMS  oechod  chat  minimizes  sample  damage  and 
allows  possible  surface  chemical  studies.  The  former  is  often  referred  to  as  dynamic  SIMS 
and  che  latter  as  static  SIMS  (43).  The  application  of  ion  beam  techniques  co  depth 
profiling  using  mass  spectromacrlc  detection  where  both  the  beam  size  and  primary  ion 
current  are  large  is  also  referred  co  as  dynamic  SIMS.  Finally,  for  the  esse  where  che  beam 
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Figure  11  -  A  general  view  of  a  secondary  ion  mass  spectrometer  system. 

From  reference  12. 

dose  is  somewhat  larger  chan  necessary  to  maintain  static  conditions,  the  term  quasi¬ 
dynamic  SIMS  has  been  applied,  tt  would  seem  that  this  approach  minimizes  the  information 
content  of  the  signal. 

We  have  learned  a  great  deal  about  the  sputtering  process  by  analyzing  the  mass,  energy 
and  ejection  angle  of  the  secondary  ions.  Of  particular  interest  is  that  in  addition  to 
ejected  monatomic  ions,  there  are  also  high  yields  of  molecular  duster  ions  chat  are  found 
to  eject.  These  were  first  observed  by  Honig  (44)  from  Ag  bombarded  by  200-400  keV  Kr+  and 
Xe+  ions  and  later  from  semiconductors  and  insulators,  graphite  (45)  and  tungsten  (46).  For 
the  case  of  tungsten,  as  shown  in  Figure  12,  clusters  containing  up  to  12  atoms  are  visible. 
Although  many  of  the  variables  which  effect  the  yields  of  cluster  ions  have  been 
investigated,  the  results  of  the  studies  have  not  provided  much  insight  into  how  the  process 
occurs.  An  exception  is  the  energy  distribution  of  the  dimer  ions  relative  to  the  monomer 
ions.  As  shown  in  Figure  13  for  niobium,  the  energy  distribution  of  the  dimers  is  more 
strongly  peaked  toward  lower  values  of  the  kinetic  energy  (47).  Similar  observations  of  the 
energy  distribution  by  Konnen  et  al.  for  ejection  of  neutral  K2  (48)  led  the  workers  to 
conclude  thet  the  dimers  were  formed  by  a  recombination  of  ejecting  monomers  above  the 
solid.  The  presence  of  molecular  cluster  ions  also  was  recognized  to  provide  potentially 
more  Information  about  the  chemlsery  of  the  solid.  Their  analysis  was  lanediecely  exploited 
in  the  ion  microprobe  (49).  The  possibility  thet  these  species  formed  after  considerable 


Figure  12  -  Mesa  spectrum  of  tungsten  clusters  sputtered  by  150  keV 
xenon  Ions .  The  energy  of  the  secondary  particle  s 
recorded  Is  30  eV,  the  energy  uldth  0.83  eV.  The  numbers 
above  the  brackets  give  the  ordinate  amplification. 

From  reference  46. 
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Figure  13  -  Secondary  Ion  spectra  of  Mb*  and  Sbj^,  not  corrected  for 
changes  In  emittance  of  the  secondary  ion  source. 


rearrangement  and  chac  Che  primary  Ion  Induced  considerable  structural  damage  in  the  sample 
loomed  as  major  factors  mitigating  against  the  dynamic  SIMS  methods  as  viable  structural 
characterization  cools. 

The  analysis  of  clusters  using  static  SIMS  at  least  removes  the  amiblgultles  Introduced 
to  the  sample  by  the  primary  beam  and  a  number  of  workers  have  now  attempted  to  exploit  this 
fact.  Sennlnghoven  (43)  first  proposed  that  for  a  number  of  different  metals  exposed  to 
oxygen,  chat  observed  clusters  of  the  type  Ma0a*  were  formed  from  contiguous  atoms  In  the 
solid  and  that  the  solid  structure  could  therefore  be  reconstructed  by  an  appropriate 
analysis  of  che  duster  Ion  Intensities.  This  Idea  was  further  reinforced  by  the  work  of 
3uhl  and  Freislnger  (30)  who  studied  ZnS  crystals  and  found  mainly  polyatomic  Ions  such  as 
ZnS*,  ZnS" ,  ZnoST*-  and  ZnSj- .  Other  polyatomic  Ions  such  as  Zn2+  or  Sj”  which  were  not 
originally  present  In  the  crystal,  ware  observed  only  In  very  low  Intensities.  Meanwhile,  It 
became  clear  that  static  SIMS  could  become  a  useful  tool  for  che  detection  of  many 
non-volatile  organic  and  Inorganic  molecules.  Sennlnghoven  first  showed  that  when  a  thin 
film  of  an  amino  acid  on  a  silver  foil  was  bombarded  with  Ar+  ions  that  a  duster  ion  equal 
to  che  molecular  weight  of  the  organic  compound  with  a  proton  attached  was  ejected  from  che 
solid  Intact  (31).  Since  chat  discovery  many  groups  have  exploited  this  observation  and  now 
such  delicate  molecules  as  proteins  with  molecular  weight  above  Z000  aou  can  be  observed. 

The  generality  of  this  assessment  has  been  questioned  both  theoretically  and 
experimentally.  For  example,  Rabalals  and  co-workers  found  that  CsCl  diluted  Into  KC1  still 
produced  significant  numbers  of  Css4,  duster  Ions  and  thac  recombination  of  ejecting 
monomers  could  occur  over  several  hundred  angstroms  (SZ).  Thus,  It  appears  thac  under 
certain  circumstances,  rearrangement  is  observed  but  in  others,  Che  dusters  are  ejected 
from  che  sample  Intact.  We  shall  see  In  lacer  sections  how  these  observations  can  be 
collected  together  to  develop  a  more  unified  theory  of  the  duster  emission  process. 

Of  further  Interest  was  che  proposal  thac  static  SIMS  could  be  utilized  to  characterize 
monolayers  formed  during  chemisorption  of  small  molecules  onto  solid  surfaces.  For  example, 
Bennlnghoven  studied  the  reaction  of  02  with  a  clean  polycryscalllne  vanadium  surface  ac 
room  temperature  (43).  As  shown  In  Figure  14,  as  che  surface  becomes  oxidized  a  number  of 
metal-oxygen  dusters  begin  to  appear.  At  high  exposure  ( 100L  O2  where  IL  -  1  Langmuir  * 
10"6  torr-sec),  the  SIMS  spectrum  Is  similar  to  that  found  for  V2O5,  Indicating  the 
formation  of  an  oxide  layer.  The  Ion  yield  fluctuations  ac  low  coverage  are  more  difficult 
to  understand,  but  It  was  felt  that  this  behavior  was  characteristic  of  a  chemisorbed  oxygen 
layer . 

Further  possibilities  for  che  study  of  monolayers  ware  put  forth  by  Barber  and 
Vtckerman  and  co-workers  who  studied  chemisorption  of  CO  on  dean  polycryscalllne  Ml 
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Figure  14  -  Chang*  of  secondary  Ion  emieslon  during  O2  exposure  of  a 
clean  vanadium  surface,  from  reference  43. 

surfaces  (53).  they  found  cluster  Ions  of  th*  type  N^4-,  N1C0+ ,  and  MliCO4,  but  none  of  the 
type  SIC4-,  HUr*-,  SljC4,  or  SlnO4-.  Indicating  that  CO  adsorbs  la  a  molecular  state  on  Ml. 
Further,  by  comparison  to  vibrational  spectra  for  this  system  they  were  able  to  make  a 
correlation  between  the  N^CO4-  Ion  and  bridge  bonded  CO,  and  the  NiCO4-  Ion  and  linear  bonded 
CO.  Their  rationale  for  this  observation  vu  that  the  molecular  cluster  was  ejected  as  a 
fragment  of  the  surface  without  rearrangement  as  schematically  shown  in  Figure  15.  These 
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Figure  15  -  Surface  processes  which  would  give  rise  to  the  appearance 
of  S1C04-  and  NlnCO41.  From  reference  53. 

workers  made  similar  conclusions  with  the  adsorption  of  sthylene  on  SI  where,  both  SliC^Ht'*’ 
and  SIC2H4*  tons  were  observed  In  the  SIMS  spectrum  (54).  Although  the  Interpretation  of 
many  of  the  abova  raaults  has  baen  questioned  In  recenc  years,  thtra  is  no  doubt  that  the 
static  STMS  approach  coupled  with  the  study  of  the  cluster  Ions  offers  a  unique  method  of 
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surface  analysis.  To  further  Che  Interpretation  of  secondary  ion  yields,  however,  ue  oust 
move  to  expand  our  cheorecical  understanding  of  the  Ion  impact  event . 

III.  Theoretical  Descriptions  of  the  Ion  3ombardment  Event 

AC  chls  polnc  In  the  discussion,  It  Is  dear  chat  many  consequences  of  the  primary  Ion 
Impact  can  be  observed  experimentally.  For  example,  ue  have  given  examples  where 
measurements  of  che  sputtering  yield  as  a  function  of  ejection  angle,  secondary  particle 
energy  and  mass  for  both  neutral  atoms  and  secondary  ions  have  been  reported.  To  fully 
appreciate  Che  significance  of  these  results  and  to  effectively  Interpret  the  fascinating 
SIMS  speccra  from  monolayers,  a  detailed  theoretical  treatment  is  clearly  required.  There 
are  essentially  two  aspects  to  chls  problem.  First,  we  must  be  able  to  determine  how  the 
nuclear  positions  of  Che  acorns  or  molecules  In  che  sample  change  with  time  In  response  to 
the  impact  of  che  primary  Ion.  Secondly,  ue  must  have  a  quantitative  description  of  the 
ionization  processes  which  cause  certain  species  to  leave  the  sample  surface  as  Ions.  In 
this  section,  we  shall  briefly  review  some  of  che  existing  approaches  to  solving  these  two 
problems  but  will  place  perhaps  undue  emphasis  on  the  models  preferred  by  the  author. 

A.  Sputtering  Field  Calculations 

There  are  basically  two  approaches  to  calculating  the  energy  dissipation  of  che 
primary  Ion.  The  first,  and  most  widely  used,  utilizes  a  statistical  model  to  solve  the 
transport  equations  for  che  momentum  deposition.  This  approach  allows  the  yield  to  be 
expressed  In  terms  of  a  simple  formula  which  can  be  easily  used  by  experimentalists.  The 
second  approach  receiving  considerable  attention  Is  to  model  the  Impact  more  explicitly 
using  classical  trajectory  methods  to  follow  the  motion  of  the  relevant  particles.  This 
method  usually  requires  a  large  computer  and  generally  does  not  provide  simple  answers. 

In  general,  there  are  a  series  of  events  which  must  be  described  following  the  Impact 
of  the  primary  Ion.  First,  Che  incoming  Ion  penetrates  the  3olld  after  striking  a  target 
atom.  This  atom  then  attains  considerable  momentum  In  the  direction  away  from  the  surface 
,nd  collides  with  neighboring  atoms  producing  secondary  displacements.  Some  of  these  recoil 
atoms  have  the  possibility  of  being  scattered  bank  coward  the  surface  with  sufficient 
kinetic  energy  to  escape  from  che  sample.  There  will  also  be  many  recoil  atoms  moving  Inside 
che  solid,  slowing  down  only  as  they  experience  more  and  more  collisions.  The  primary  Ion 
and  the  target  atom  may  either  come  to  rest  inside  the  solid  (implantation)  or  be  ejected 
themselves.  Thus,  the  sputtering  process  is  a  complex  event  involving  multiple  atomic 
colllsona  (or  cascade)  initiated  by  che  primary  ion. 

1 .  Tranaport  Theories 

There  have  been  a  number  of  attempts  to  derive  analytical  expressions  which  predict 
various  properties  of  the  sputtering  phenomenon.  The  most  widely  used  are  those  of  Slgmmd 
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(55-57)  whose  yield  formula  are  ubiquitous  throughout  the  sputtering  literature  and  of 
Thompson  (58)  who  has  provided  a  relationship  for  the  yield  as  a  funccion  of  secondary 
particle  velocity.  In  both  of  these  approaches,  the  Initial  assumption  is  that  each  target 
atom  generates  one  secondary  collision  cascade  In  which  energy  is  shared  by  a  series  of 
collisions.  Next,  It  Is  assusied  that  the  collisions  themselves  occur  only  between  cwo  atoms 
at  a  time.  These  approximations  have  a  number  of  important  consequences.  First,  any 
electronic  energy  loss  processes  are  ignored.  This  is  probably  a  realistic  approximation 
for  metals  bombarded  at  less  chan  50  kaV  by  heavy  Ions,  although  for  insulators  such  as 
alkali  halides  (59),  coulomblc  interactions  may  contribute  to  sputtering.  For  elastic 
collisions,  the  assumption  chat  each  collison  sequence  can  be  created  separately  restricts 
application  of  the  theory  to  the  region  of  linear  collision  cascades  (10,55).  Aa  shown  in 
Figure  16a,  this  condition  applies  only  when  a  small  fraction  of  all  atoms  within  a  certain 
volume  are  in  motion.  When  more  than  one  cascade  Is  Initiated  by  the  primary  Ion,  as  shown 
in  Figure  16b,  or  when  a  large  fraction  of  Che  atoms  near  the  Impact  point  are  moving,  the 
approach  begins  Co  breakdown.  This  region  is  often  referred  to  as  a  dense  cascade  or  a 
spike.  One  final  consequence  is  that  the  binary  collision  approximation  (BCA)  should  only 
apply  to  collisions  which  exceed  several  hundred  eV.  For  example.  In  Figure  16a,  the  sizes 
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FljutaJJs  -  (a)  Linear  collision  cascade.  The  structure  is  preserved,  and 

a  small  fraction  of  the  stoma  Is  In  motion. 

(b)  Dense  cascade  (spike).  The  structure  Is  destroyed 
locally.  All  acorns  within  the  spika  volume  are  In  motion. 

From  reference  10. 
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of  che  acomt  are  drawn  such  chac  lc  Is  Impossible  for  more  chan  two  acoae  Co  be  lnceracciog 
simultaneously .  If  che  atomic  radii  had  been  drawn  such  chac  adjacent  atoms  nearly  touched, 
as  in  che  case  for  collision  energies  below  —10  aV,  >en  binary  collisions  are  Indeed 
improbable.  It  is  disturbing  that  moat  of  the  collisions  chat  lead  to  the  ejection  of  atoms 
in  che  solid  are  chose  whose  energy  is  below  10  eV  as  illustrated  by  che  secondary  particle 
energy  distributions  given  in  Figures  8  and  10. 

In  che  Sigmund  theory,  che  solid  is  approximated  by  an  amorphous  array  of  atoms  chac 
inceract  wlch  differencial  cross  sections  developed  for  calculations  aimed  coward  predicting 
the  range  of  energetic  atoms  in  solids  CSS).  The  nuclear  stopping  cross  section  Sn(E),  can 
be  evaluated  from  these  quantities .  The  flow  of  energy  through  the  cascade  is  Chen 
determined  using  che  Boltzmann  transport  equation.  The  sputtering  yield  is  evaluated  from 
the  flux  of  particles  crossing  an  imaginary  plane  at  a  point  in  che  solid  where  che  target 
acorn  is  given  its  initial  energy.  To  Include  che  effect  of  che  surface,  lc  is  assumed  chat 
che  arom  must  surmount  a  potential  barrier,  U,  acting  on  the  perpendicular  component  of  the 
velocity.  This  potential  will  act  to  keep  slow  moving  particles  from  escaping  the  barrier 
and  will  deflect  the  faster  moving  ejecting  particles  away  from  the  surface  normal.  This 
correction  is  termed  che  planar  surface  binding  model  (55,56).  If  che  collision  cascade  is 
restricted  to  the  linear  regime  by  utilizing  bombardment  energies  below  a  few  keV  and  by 
utilizing  light  ion  and  target  masses,  the  equation  for  Che  sputtering  yield  becomes 

5  .  C  SqiE)  (g) 

U 

where  C  is  a  constant  for  a  given  primary  ion  and  energy  and  target  material.  The  simpli¬ 
city  of  this  result  -  chat  the  sputtering  yield  is  proportional  to  the  nuclear  stopping 
power  of  the  solid  and  inversely  proportional  to  the  surface  binding  energy  -  accounts,  in 
part,  for  the  popularity  of  the  Sigmund  theory.  In  fact,  it  hae  been  rather  successfully 
applied  in  a  number  of  Instances  as  has  been  reviewed  (60).  The  direct  correlation  between 
che  stopping  power  and  the  yield  is  particularly  striking  as  is  illustrated  in  Figure  17  for 
Ar+  on  Cu  (56,61).  The  relationship  between  yield  and  surface  binding  energy  (or 
sublimation  energy  of  the  solid)  has  been  ueed  to  explain  differential  sputtering  (62)  and 
variations  in  sputtering  yield  acroes  the  periodic  cable  (55,63).  The  cheory  has  also  been 
applied  co  situations  where  lc  should  noc  work,  probably  since  it  has  been  che  only  reel 
theoretical  approach  available  co  experimentalists . 

A  similar  theoretical  modal  has  been  developed  by  Thompson  to  explain  the  energy 
distribution  of  ejected  atoms  (58).  His  model  develops  an  expression  for  che  flux  of  scorns 
Inside  che  solid  crossing  any  internal  surface  plane.  For  an  ejecting  particle  of  kinetic 
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Figure  17  ~  Sputtering  yield  as  a  function  of  the  primary  ion  energy.  The 
lines  are  calculated  yields  from  Sn(E)  while  the  points  are 
experiment  curves  taken  from  various  groups.  From  references 
55  and  56. 

energy  E,  this  flux  is  given  as 

XE)dEdJ>(C'Ecos8/(E+C)n+l]  '(2frsiud)d8dE  (7) 

where  XE)  is  the  flux  density  energy  distribution  with  units  of  number  per  unit  ares  per 

unit  time  per  unit  energy.  The  solid  angle  d!)  «  sin8d6dp  »  2*sln6d9;  dp  is  the  azimuthal 

angle  about  which  the  flux  is  assumed  to  be  isotropic.  According  to  Thompson  (58),  n~2  and 

C  Is  a  constanc  for  a  given  primary  ion,  energy  and  target  material.  The  appearance  of  ’J 

in  Equation  7  is  a  consequence  of  the  incorporation  of  the  planar  surface  binding  model  in 

the  cheoreclcal  development.  For  E»0,  the  energy  distribution  should  call  off  as  E~2.  The 

most  probable  energy  should  occur  at  approximately  U/2.  There  have  been  many  experiments 

aimed  coward  casting  this  equation.  Probably  the  most  reliable  data  so  far  are  produced 

from  the  laser  Induced  fluorescence  experiments  on  sputtered  uranium  and  zirconium  as  shown 

in  Figure  10  (64,65).  In  this  case,  the  value  of  a  in  Equation  7  was  found  to  be  between 

1.75  and  2.0  with  U  *  5.2  eV.  The  development  of  the  Thompson  theory  waa  a  major  advance  in 

our  understanding  of  sputtering  and  stood  to  elucidate  the  differences  between  sputtering 

and  simple  evaporation.  In  the  latter  case,  the  energy  distribution  should  appear  to  b« 
Maxvell-Eoltzaan  like  in  shape. 


Although  boch  of  these  theories  were  helpful,  they  were  not  well-suited  to  studying 
crystalline  solids  and  angular  anisotropies.  They  were  forced  to  utilize  the  3CA  and  to 
employ  oversimplified  Interaction  potentials  to  obtain  analytical  results.  It  was  not 
possible  to  explain  the  factors  that  Influence  cluster  formation  except  by  a  clever  analysis 
of  the  energy  distribution  formula  (66-68).  The  theories  have  been  no  help  In  understanding 
the  basic  phenomena  In  SIMS. 

2 .  Molecular  Dynamics  Calculations 

Many  of  thesa  difficulties  can  be  overcome  by  utilizing  a  molecular  dynamics 
calculation  on  a  large  ensemble  of  atoms  to  compute  actual  nuclear  positions  as  they  change 
In  dme  subsequent  to  the  primary  Ion  event.  Classical  dynamics  calculations  have,  of 
course,  been  very  successful  In  explaining  trajactorles  In  acom-diatom  scattering  (69), 
properties  of  liquids  (70)  and  even  the  solvation  of  large  molecules  like  dlpeptldes  (71). 
For  describing  the  sputtering  process  this  approach  has  the  distinct  advantage  of  utilizing 
many  fewer  approximations  chan  required  for  the  statistical  theories.  On  the  the  other 
hand,  no  simple  equation  falls  out  of  Che  calculations,  although  Important  concepts  ma^r 
emerge  from  the  resulting  numbers.  The  calculations  often  prove  very  useful  in  testing  the 
validity  of  possible  analytical  theories.  In  general,  these  calculations  require 
considerable  computer  time.  With  the  recent  surge  In  computer  efficiency.  It  would  appear 
this  Is  becoming  leas  of  a  difficulty. 

The  first  efforts  at  developing  the  molecular  dynamics  procedures  were  undertaken  by 
Vineyard  and  co-workers  during  the  early  1960's  (72-74).  These  workers  constructed  a  model 
microcrystalllce  consisting  of  several  hundred  atoms.  A  two-body  central  Interatomic 
potential  was  used  to  describe  the  interaction  between  each  pair  of  atoms  as  energy  was 
dumped  into  Che  system.  The  lengthy  time  required  for  the  computation  restricted  the 
application  to  mechanistic  studies  of  how  atoms  moved  in  the  solid  and  could  noc  be  compared 
directly  to  experiment.  Subsequently,  Harrison  (75,76)  was  able  to  perform  a  number  of  more 
lengthy  calculations  for  Ar*"  Ions  bombarding  a  Cu  single  crystal.  Vlelds  and  angular 
distributions  could  be  calculated  chat  were  not  coo  far  away  from  the  experimental  results. 
Hsrrlson  was  conclusively  able  to  show,  for  txsmple,  thee  Wehner’s  spot  pectems  were  due 
mainly  to  structural  effseca  of  the  top  etomlc  layer  of  the  metal  (77),  rather  than  from 
focussed  collisions  from  deep  within  the  solid  ss  had  been  previoualy  believed  (78). 

Similar  mechanistic  studies  to  dectrmlne  the  types  of  scattering  events  observed  in  metals 
were  also  undertaken  by  Jackaon  (79,30). 

To  shorten  the  required  computer  time,  Robinson  and  co-workers  developed  a  model  that 
Incorporated  the  BCA,  the  seme  approximation  as  that  Implemented  In  the  energy  transport 
theory  (31,82).  With  this  approximation,  these  workers  were  able  to  follow  a  series  of 
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collisions  within  Che  cascade  by  only  calculating  che  energy  transferred  during  each 
collision.  The  cooputer  tine  was  reduced  by  many  orders  of  magnitude  which  allowed  direct 
coaparlson  to  a  number  of  experimental  yield  measurements.  Although  this  type  of 
calculation  is  well-suited  to  modeling  IcaV  ion  scattering  from  solids  where  the  BCA  is 
applicable  (83),  there  are  serious  questions  as  we  have  mentioned  earlier,  about  che 
validity  of  using  it  in  sputtering  calculations.  Boblnson  has  recently  recognized  this 
difficulty  and  incorporated  an  adjustable  parameter  into  his  computer  code  which  accounts 
for  multiple  collisions  (81).  The  latest  attempts  to  compare  yields  and  energy 
distributions  to  experimental  values  of  Ar+  ions  on  amorphous  Cu  have  been  reasonably 
successful  (85).  In  this  work,  che  authors  also  Incorporated  both  che  planar  surface 
binding  model  and  isotropic  surface  binding  model  into  their  calculation.  In  the  latter 
case,  che  ejecting  particle  losses  energy  a  regardless  of  its  ejecting  direction.  It  is  not 
clear  yet  which,  if  either,  of  these  two  approaches  to  describing  the  atomic  trajectory  as 
it  crosses  che  surface  is  best.  The  calculated  yields  are  higher  by  about  30X  using 
isotropic  binding,  but  since  the  value  chosen  for  0  is  somewhat  arbitrary,  it  would  appear 
this  is  a  meaningless  difference.  The  angular  distributions  calculated  for  the  two  surface 
binding  models  should  be  quite  different  since  a  bending  of  che  ejecting  atom  away  from  the 
surface  normal  occurs  only  for  che  case  of  che  planar  binding  model.  Unfortunately, 
detailed  comparisons  to  experiment  have  not  yet  been  successful  in  distinguishing  between 
che  two  approaches.  Other  groups  have  developed  similar  models  using  the  BCA  (86-89), 
although  none  has  been  successful  in  explaining  most  sputtering  experiments. 

The  advent  of  che  present  day  computer  has  allowed  vineyard's  original  ideas  to  be 
exploited  to  a  much  greater  extent.  During  the  remainder  of  this  chapter  we  shall  focus  on 
these  calculations  since  they  have  been  remarkably  successful  at  predicting  a  wide  range  of 
experimental  observables  found  tor  the  sputtering  of  neutral  atoms.  In  addition,  there  are 
some  real  possibilities  for  coupling  these  types  of  calculations  which  describe  the  nuclear 
motion  of  the  atoms  to  chose  which  describe  che  electronic  motion  as  we  shall  see  shortly. 

The  computation  of  classical  trajectories  using  the  molecular  dynamics  procedure  rests 
on  Newton's  equaclon  of  motion.  For  a  particle  1  of  mass  m^,  the  equation  of  motion  is 

mj?i  -  ?•  (8) 

^  Ia  ^ 

where  r^  is  the  position  vector,  r^  is  the  acceleration  and  F^  the  force.  For  N  particles 
there  will  be  3N  coupled  second  order  differential  equations  chat  need  to  be  solved. 
Computationally,  it  is  more  convenient  to  utilize  Hamilton's  form  of  che  equations  of  motion 

as 

m^vj  *  (9) 


rl  ’  V1 


(10) 
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where  is  the  particle's  velocity,  since  there  are  now  6N  coupled  first  order  differential 
equations.  These  equations  provide  the  position  and  the  time  dependence  of  the  velocity  and 
are  easier  to  solve  than  the  second  order  equations,  For  an  arbitrary  value  of  F*i  the  atonic 
moclon  can  be  determined  using  a  finite  difference  algorithn  and  an  lncrenental  tine  step  At 
(90) .  For  example,  the  average  force  between  c  and  t+At  Is 

<?!>  -  l/2[F'1(fi*(t+dt)  +  F>iC?l<t))!  (11) 

where  Is  the  predicted  position  at  t+At  and  is  given  by 

^‘(t-t-At)  -  ri.(t)  +Vi(t)At  +  Fifrtlt)  )At2/2ni  (12) 

Once  <Fi>  Is  known,  the  corrected  positions  and  velocities  at  t+At  are  then 

1?l(t+At)  »  'rj(t)  +  ?i(t)At  +  <Fi>At2/2ni  (13) 

and  v^t+At)  -  vi(t)  +  <Fi>At/mi  (14) 

For  a  single  component  system  the  value  of  At  Is  computed  from  the  value  of  the  moat 
energetic  particle  remaining  in  the  solid  at  time  t.  For  particles  of  disparate  mass,  this 
Is  not  satisfactory  since  lighter  particles  move  farther  than  heavier  ones  with  the  same 
kinetic  energy.  A  better  scheme  for  selecting  At  Is  to  base  the  determination  on  the 
largest  difference  between  the  predicted  and  corrected  position  Ar.  In  this  case  (91) 

&tnew  *  Atol(l(e/Ar)a  (15) 

where  c  Is  an  Input  parameter  of  order  10~3  to  10~4  chosen  to  yield  the  desired  Integration 
accuracy  of  the  final  observables.  The  value  of  a  has  been  arbitrarily  set  to  0.3  for 
similar  reasons  (92).  This  approach  allovs  the  tlmestep  to  vary  during  calculation  of  the 
■trajectory.  At  the  beginning  of  the  trajectory  the  tlmestep  Is  typically  ~10-16  s  and 
eventually  Increases  to  -10"14  s.  A  tocal  of  100-200  Integration  steps  are  required.  It  Is 
of  interest  that  for  these  systems  a  low  order  predictor-corrector  Is  most  efficient 
computationally  since  the  forces  vary  rapidly  with  distance. 

In  all  calculations  to  date,  the  Initial  coordinates  of  the  atoms  are  set  at  their 
equilibrium  position  with  zero  Initial  velocity.  As  shown  In  Figure  13,  the  primary  Ion  Is 
given  a  specific  kinetic  energy  E,  with  polar  angle  9  and  azimuchal  angle  a  such  that 


v  »  (2E/m)^2  ■  (vx*  ♦  Vy2  +  Vj2)2^2  (16) 

vt  •  -v  cos  9  (17) 

vx  •  v  sin  9  cos  t  (18) 

Vy  •  v  sin  9  Sin  >  (19) 
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Figure  18  -  Coordinate  system  used  In  determining  angular  distributions. 

The  force  la  calculated  from  the  gradient  of  the  potential  energy,  V, 

as  -  -V^V  <ri,?2...r*n>  (2°) 

For  a  nany-body  solid,  we  would  need  a  aulti-dlmenslonal  potential  surface  that  accounted 
for  all  possible  positions  of  che  particles.  With  our  presenc  knowledge  of  the  solid  state, 
this  potential  surface  Is  unknown  buc  Is  usually  approximated  using  a  sum  of  palr-wlse 
additive  potentials  as 

7  •  I  I  7lj  (Uj)  (21) 

l-l  j>l 

where  r^j  Is  the  distance  between  particle  1  and  j . 

The  subsequent  atomic  motion  due  to  the  primary  Ion  Impact  Is  computed  for  enough  time 
such  that  no  particles  have  sufficient  kinetic  energy  to  escape  che  solid.  For  a  metal  such 
as  copper,  we  have  found  that  after  the  most  energetic  particle  has  less  chan  2  ev  of 
kinetic  energy  chat  Integration  for  longer  times  does  not  produce  any  more  ejecced 
particles. 

After  termination  of  the  calculation,  the  particles  thac  have  ejected  are  tested  for 
possible  multlmer  formation  (93,91).  To  check  for  che  formation  of  sulclmers,  one  merely 
computes  the  relative  kinetic  energy  Zr,  plus  the  potential  energy  7ij ,  for  all  pairs  of 
ejected  atoms  1  and  j .  If  che  tocal  energy  of  the  dimer 

dimer  dimer 

Scot  -  =r  +  7lj  (22) 

Is  negative,  then  the  tested  dimer  Is  considered  to  be  bound.  Further,  It  Is  often  found 


that  several  bound  dimers  are  formed  above  the  surface.  If  some  of  these  have  common  atoms', 

l.e.  If  the  dimers  are  linked  or  overlapping,  then  higher  multlaers  may  possibly  exist. 

cluster 

If  this  condition  Is  found,  Etot  where 


cluster  “pi 


♦  T  \ 


tot  r  1-i  j<i  J 

with  a  being  che  number  of  atoms  In  the  cluster.  Is  recalculated  for  all  the  atoms  In  the 

cluster 

linkage.  As  In  che  dimer  analysis.  If  Etot  for  the  atoms  In  the  linkages  la  less  than 

zero  then  che  atoms  are  considered  to  be  a  duster. 

There  are  other  possible  definitions  of  what  constitutes  a  cluster.  The  requirement 

thac  each  atom  In  the  mulcimer  be  bound  Co  another  atom  may  be  an  overly  stringent 

duster 

requirement  for  cluster  stability.  Rigorously,  any  collection  a f  n  atoms  with  &toc  <0 
Is  considered  temporarily  bound.  Since,  as  ve  shall  see,  there  are  uncertainties  in 
It  is  prudent  to  use  che  conservative  definition  of  a  cluster.  Further,  for  i>2  It  Is 
possible  chat  the  cluster  will  decompose  before  reaching  che  detector.  For  che  higher 
dusters,  then,  che  model  really  only  tests  to  see  If  there  are  significant  numbers  of 
ejecting  atoms  that  are  in  a  favorable  spaclal  position  with  low  enough  relative  kinetic 
energy  to  experience  binding  interactions. 

In  addition  co  evaluating  the  possibility  of  duster  formation,  che  classical  dynamics 
procedure  allows  all  possible  observable  properties  of  the  ejected  neutral  atoms  to  be 
calculated.  The  yield  Is  determined  by  counting  the  number  of  particles  which  are  moving 
away  from  the  surface  and  are  not  bound  to  Che  solid.  The  kinetic  energy  of  the  ejected 
atoms  (or  dusters)  Is  computed  from  the  known  velocities  as 


E  •  l/Zm  (vx^  +  Vy2  +  Vj2) 


The  angular  distributions  are,  of  course,  known  from  the  velocity  components  since 


cos  9  »  -vz/v 

(25) 

tan  $  «  Vy/vx 

(26) 

It  is  also  possible  co  trace  che  path  of  Individual  atoms  as  they  move  through  the  cascade 
volume.  As  we  shall  see,  this  tedious  job  is  of  extreme  interest  in  evaluating  the 
Importance  of  proposed  ejection  mechanisms.  For  example.  It  Is  possible  to  say  with 
certainty  whether  multiple  collisions  dominate  the  sputtering  process  or  whether  there 
exists  a  sec  of  experimental  conditions  where  che  linear  cascades  (Figure  16a)  are  Indeed 
observed  to  occur. 


The  major  difficulty  with  this  model,  or  any  modal,  of  the  loo/ solid  Interaction  Is 
chat  urn  have  only  a  Halted  knowledge  of  the  forces  acting  between  the  aeons.  Sons 

Information  Is  available.  For  example,  gas  phase  scattering  axperlaenta  at  a  variety  of  _ 

energies  can  be  used  to  aap  out  Che  Interaction  potential  surface  (95).  Extrapolating  this 
Information  to  the  solid  state  la  probably  reasonable  above  energies  of  a  few  hundred  eV 
since  the  solid  atoms  are  essentially  not  Interacting  anyway  at  these  energies.  Other 
approaches  have  Included  attempts  to  fit  various  forms  of  a  potential  to  observables  such  as 
elasticity  and  heat  of  subllaatlon  (95,96). 

Whenever  che  BCA  is  not  used,  pair-vise  additive  potentials  have  always  been  assumed. 
This  approximation  needs  to  be  Improved  for  two  reasons.  First,  many  solids  have  directional 
bonds  -  Si,  CaAs  or  ice  are  examples  where  this  correction  needs  to  be  taken  into 

account.  Second,  there  Is  a  major  problem  in  attempting  to  move  an  atom  from  the  bulk  of 

che  solid  where  It  may  have  12  neareat-neighbors  to  vacuum  where  it  exists  as  an  Isolated 

atom.  In  the  calculations,  this  difficulty  has  been  Ignored  except  when  analyzing  for 

duscers  where  che  potential  used  In  the  solid  is  switched  to  a  potential  more  appropriate 
for  a  gas  phase  species.  A  rigorous  description  of  che  actual  ejection  process,  then,  is 
lacking  in  all  models  of  sputtering  yet  proposed. 

For  a  single  component  metal  surface  such  as  Cu,  che  Cu-Cu  Interaction  is  often 
represented  using  the  interaction  potential  suggested  by  Harrison  (76).  This  potential 
utilizes  a  3orn-Hayer  exponential  function  to  approximate  che  high  energy  collisions  (>20 
eV)  and  a  Morse  potential  with  parameters  determined  by  Anderman  (96)  to  approximate  the 
low  energy  (<5  eV)  collisions.  A  cubic  spline  is  used  to  connect  che  two  parts.  The  precise 
form  of  the  interaction  between  atoms  1  and  j  separated  by  distance  r^j  is 


7ij 

■  Ae  “®rij 

rij<ra 

(27) 

7U 

*  Co  +  Cl  rij  +•  C2  rtj2  s -  C3  Tij3 

ra<rlj  <rb 

(28) 

vu 

•  De  exp(-3(rtj-re)  1[exp{-S(r2j-ra)  1-2  ] 

rb<rlj<rc 

(29) 

vu 

-  0 

rij>rc 

(30) 

where  De  is  che  well-depth  of  the  potential  and  ra  is  che  equilibrium  separation  distance. 
The  Incorporation  of  an  attractive  Interaction  helps  to  alleviate  a  number  of  the 
ambiguities  of  describing  the  atom  trajectory  as  it  crosses  che  surface.  For  example,  in 
che  "Isotropic"  and  "planar"  surface  binding  model,  a  purely  repulsive  potential  is  used 
even  during  the  ejection  process.  Then,  the  kinetic  energy  needed  to  overcome  the  surface 
binding  is  subtracted  at  the  end  of  che  calculation.  In  che  full  dynamics  calculations 
these  interactions  are  calculated  during  che  ejection  process  obviating  the  need  for  either 
the  planar  or  isotropic  binding  models. 


There  are  several  approaches  co  describe  che  Interaction  between  the  primary  Ion  and 


the  solid.  An  accurate  description  of  this  Interaction  Is  Important  only  at  fairly  high 
energies  since  moat  of  the  momentum  is  leparted  to  the  solid  after  only  a  few  collisions. 
?robably  the  aost  popular  fora  of  this  Interaction  is  the  Mr’.iere  potential  (85) 

3 

V1J  ”  zlZ2«2/rljn3j  “a**7  (-Bnrlj/*>  (31) 


where  Zje,  Z2e  are  the  nuclear  charges  of  the  projectile  and  target  atoms,  a  Is  a  screening 
length,  a  •  (0.35,  0.55,  0.10)  and  B  •  (0.3,  1.2,  6.0).  This  sua  at  three  exponentials  la  a 
fit  to  che  Thomaa-Fermi  potential.  Values  of  the  screening  length -are  calculated  from  the 
Firsov  formula 


a  •  (9vz/ 128) 1/3  a0(Zj1/2  ♦  Z21/2>-2/3  (32) 

where  a0  Is  Che  Bohr  radius.  It  has  been  proposed  chat  a  value  of  a  which  Is  20Z  of  Che 
value  calculated  In  Equation  32  gives  slightly  better  agreement  with  experiment.  In  his 
early  work,  Harrison  (76)  employed  the  Bom-Mayer  fora  of  the  lon-mecal  interaction  as 


7lj  ■  A  «'Brij 

rlj<re 

(33) 

7lj  -  o 

rij>ra 

(34) 

This  form  apparently  works  well  below  Ar+  energies  of  ~1  KeV  but  becomes  much  too  large  at 
energies  above  this  value.  Two  forms  of  the  Interaction  potential  for  che  Ar+-Nl  system  are 
plotted  In  Figure  19. 

It  Is  Important  to  be  sure  that  there  are  a  sufficient  number  of  atoms  In  the 
olcrocrystalllte  such  that  che  calculated  observable  quantity  does  not  change  as  the 
alcrocrystalllte  size  is  Increased  further.  In  practice,  considerable  testing  with  various 
size  models  Is  important  to  verify  this  condition.  For  600  eV  Ar*  Ion  bombarding  Cu(00l)  at 
normal  Incidence  approximately  &  layers  of  atoms  with  60  acoms/layer  are  required  for 
convergence  of  most  of  the  desired  quantities.  As  che  kinetic  energy  of  the  Ar+  Ion  Is 
Increased  further,  however,  che  number  of  acorns  thac  are  needed  Increases  proportlonat  ely 
such  that  approximately  2000  atoms  must  be  Included  to  obtain  convergence  at  5  keV  Ar+  Ion 
energies.  The  calculated  yields  for  soma  different  sized  crystals  at  1000  and  2000  eV  are 
shown  In  Table  I.  Slightly  larger  crystals  ars  also  required  for  oblique  Incident  angles. 
Much  smaller  cryetals  may  be  employed  If  only  the  higher  energy  ejected  particles  are  of 
Interest  (97). 


Table  I 

Effect  of  Crystal  Size  on  Calculated  Yield  for 
1000  and  2000  eV  Ax+  Ion  Bombardment  of  Cu(001) 

1000  eV 


Crystal  Size 
t  of  Atoms/Layer 

#  Layers 

Relative  Tield(a> 

84 

4 

0.84 

112 

4 

0.96 

84 

5 

0.90 

112 

5 

0.95 

84 

6 

0.93 

112 

6 

1.00 

144 

6 

0.92 

2000  eV 

84 

5 

0.73 

112 

5 

0.81 

144 

6 

0.97 

180 

6 

1.00 

(•^Obtained  using  a  Mollere  potential  for  the  Ar+-Cu  interaction 
and  a  Born-hayer  and  Morse  potential  for  the  Cu-Cu  interaction. 
Yields  are  normalized  to  Che  maximum  value  in  each  sec. 


In  this  calculation,  the  only  mechanism  by  which  the  crystal  can  lose  energy  is  by 


ejection  of  particles  either  from  the  surface  or  at  the  sides.  In  Harrison's  original  model, 
considerable  testing  showed  that  it  was  most  desirable  chat  any  particle  ejecced  from  the 


sida  or  che  bottom  be  allowed  to  carry  away  Its  momentum,  preventing  It  from  further 
contributing  to  the  action.  If  boundary  conditions  are  eaployed  which  artlflcally  Weep  the 
energy  within  the  crystal,  the  only  aechanlsa  for  energy  lose  Is  by  the  ejection  of  aaterlal 
through  che  surface.  Slailar  difficulties  exist  If  periodic  boundary  conditions  are 
eaployed.  The  'edge*  effects,  Chen,  are  controlled  by  asking  the  alcrocrystalllte  large 
enough  so  that  its  size  does  noc  Influence  the  observable  and  by  allowing  energetic  atoos 
moving  through  che  sides  or  che  boctoa  co  simply  carry  away  their  excess  energy. 

The  details  of  che  trajectory  are  strongly  Influenced  by  the  point  at  which  che 
primary  projectile  strikes  che  target  aeon.  For  example.  If  che  primary  ion  impacts  che 
target  atoa  head  on.  It  Is  driven  loco  che  solid  and  no  particle  ejection  Is  observed.  For 
certain  Impact  points  away  from  this  condition,  however,  tremendous  motion  Is  observed 
throughout  the  crystallite  producing  a  large  number  of  ejected  atoms.  This  distribution, 
sometimes  termed  Che  number  of  atoms  sputtered  per  incident  Ion  (ASI)  (98)  for  600  eV  Ar+  on 
Cu(OOl)  can  vary  from  zero  co  greater  than  20. 

The  fact  thac  the  yield  Is  a  sensitive  function  of  che  impact  point  coordinates 
requires  sampling  of  all  symmetrically  lnequlvalent  surface  sites .  In  Figure  20 ,  for 
example,  are  shown  parts  of  the  (001),  (110)  and  (111)  surfaces  of  an  fee  metal  such  as  Cu. 
The  zone  of  Irreducible  symmetry  Is  illustrated  for  each  case.  In  practice  somewhere 
between  100  and  1000  trajectories  are  computed  at  Impact  points  over  this  zone  to  obtain  the 
macroscopic  yield  of  particles  and  ocher  observables.  The  higher  number  of  trajectories  are 
required  If  che  yield  at  a  particular  angle  or  energy  is  required.  Detailed  analysis  of 


Figure  20  -  Crystal  surface  and  Impact  zone  for  the  three  low  index  faces  of  an  fee  mecal. 

ejected  duscers  also  demands  a  higher  number  of  trajectories  since  their  yield  Is  usually 
only  a  few  percent  of  thac  of  che  monomer.  The  statistical  reliability  of  any  result  may  be 
estimated  as  n0ps  t  2  (ngi,,)1^  where  nobt  Is  the  total  number  of  observed  evencs  summed 
over  all  the  computed  Impact  points. 
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The  dynamics  of  more  complex  model  mlcrocryitelllcet  1*  also  feasible  Co  study  using 
che  compucer  model.  For  example,  chemisorbed  scorns  on  mecals  (99),  chemisorbed  molecules 
(100,101)  on  mecals,  allcall  halides  (102)  and  alloys  (103)  have  all  been  Investigated  with 
varying  degrees  of  detail.  For  the  case  of  atomic  adsorbaces,  che  effect  of  chelr  coverage 
and  bonding  geometry  can  be  examined  to  determine  ejection  mechanisms  and  che  usual 
experimentally  measurable  quantities.  The  placement  of  atosts  on  a  clean  surface,  however, 
alters  che  Impact  tone  of  Irreducible  symmecry.  For  example,  a  single  atom  may  be  placed  In 
a  linear  bonded  or  atop  geometry  on  che  center  substrate  atom.  This  configuration 
approximates  the  zero  coverage  limit  in  that  Che  adsorbate-adsorbate  scattering  mechanisms 
cannoc  occur.  The  Impacc  zone,  however,  grows  to  1/8  of  the  entire  surface.  Other 
configurations  such  as  che  p(2x2)  (252  of  a  monolayer)  and  the  c(2x2)  (502  of  a  monolayer) 
ordered  overlayer  structures  for  three  different  adsorption  site  geometries  are  shown  in 
Figure  21.  Ho  adsorbate  atoms  have  been  placed  on  che  edge  of  the  microcrystallite  since 
these  positions  exhibit  unrealistically  wait  binding  forces.  The  Impact  zones  for  all  of 
che  cases  are  also  shown  In  Figure  21;  note  as  the  symmetry  decreases,  the  size  of  the 
impact  zone  Increases. 

The  Interaction  potentials  are  formulated  In  much  the  same  manner  as  for  the  clean 
metal  case.  If  atomic  locations  are  known  from  low  energy  electron  diffraction  or  other 
methods  sensitive  to  surface  structure  than  che  vaxue  of  re  In  Equation  29  is  known.  Also, 
the  value  of  De  can  be  estimated  from  the  heat  of  adsorption.  The  parameters  used  to 
construct  the  repulsive  potential  wall  should  only  be  a  function  of  the  atomic  number  and 
not  of  che  electronic  nature  of  the  bonding  since  che  energies  are  much  larger  chan  typical 
bond  strengths.  Many  of  che  predictions  of  che  classical  dynamics  calculations  for  these 
systems  will  be  discussed  in  future  secclons. 

B.  Ionization  of  the  Ejected  Atoms 

There  is  clearly  a  greac  deal  of  Information  concerning  che  ejection  of  neutral  atoms 
from  solids.  The  cheory  Is  in  reasonably  good  shape,  especially  since  many  of  the 
observables  can  be  quantitatively  predicted  using  the  molecular  dynamics  calculations.  All 
chat  Is  needed  is  a  reasonably  reliable  Interaction  potential  which  can  usually  be 
constructed  from  fairly  simply  determined  physical  properties.  The  question  of  lonlzadon 
at  surfaces  Is  an  entirely  different  matter.  For  this  slcuatlon,  there  have  been  literally 
dozens  of  models  proposed  with  agreement  to  experiment  usually  no  beccer  than  a  couple  of 
orders  of  magnitude.  There  has  been,  therefore,  a  focus  on  developing  rather  empirical 
experiments  and  correlations  to  try  co  get  this  aspect  of  the  sputtering  phenomenon  under 


control 
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Figure  21  -  Possible  placements  of  the  oxygen  atoms  on  the  microcrystallite. 

The  shaded  circles  represent  oxygen  atoms  and  the  open  circles 
represent  substrate  atoms.  The  Impact  zones  are  also  shown. 

To  provide  some  justification  for  the  development  of  the  exisclng  theories  which 
actempt  to  explain  surface  ionization,  it  is  useful  to  lay  ouC  the  important  differences 
between  atom  yields  and  ion  yields  from  the  same  substrace.  In  the  moat  general  sense,  the 
neutral  yields  vary  only  by  as  much  as  a  factor  of  10  for  a  given  primary  Ion  energy  for 
virtually  all  elements  of  the  periodic  table.  As  shown  In  Figure  22  for  both  positive  and 
negative  Ions,  however,  It  is  seen  that  may  vary  by  over  5  orders  of  magnitude  (104). 
Next,  there  appears  to  be  a  rough  Inverse  correlation  between  Ionization  efficiency  and 
increasing  lonlxation  potential  for  positive  ions  and  increasing  electron  affinity  for 
negative  Ions  (105).  The  neutral  yields  exhibit  no  such  correlation.  And  finally,  there 
has  been  a  strong  dependency  noted  between  the  electronic  properties  of  the  surface  (l.e. 
work  function)  and  che  magnitude  of  R+  for  a  given  element.  For  example,  the  yield  of  Ni+ 
(as  well  as  Nlj*  and  Nl]"*-)  from  a  clean  Ni  surface  bombarded  by  1  keV  Ar+  is  near  zero 
(106),  whareas  the  yield  of  Nl  atoms  Is  approximately  4(1),  Upon  formation  of  50?  of  a 
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Figure  22  -  (a)  Relative  secondary  positive  ion  yield  (M*)  due  to  0" 

bombardment  vs.  atomic  number  of  secondary  ion.  (e)  Pure 
element,  (&)  Compound. 

(b)  Relative  secondary  negative  ion  yield  (M~)  from  Cs* 
bombardment  vs.  atomic  number  of  secondary  ion.  Relative 
intensity  is  given  in  units  counts/ second/ 10“9  A.  (e)  Pure 
element.  (4)  Compound.  B.D.  •  Barely  Detectable.  S.D.  - 
Hot  Detected.  From  reference  105 . 


monolayer  of  oxygen,  hovever,  the  Hi+  yield  increases  by  several  orders  of  magnitude  such 
that  K*'~0 .1  (106).  The  sensitivity  of  the  Intensity  of  the  Ni+  signal  to  oxygen  coverage  is 
illustrated  in  Figure  23.  It  is  this  tremendous  variation  which  causes  so  much  difficulty 
in  attempting  to  obtain  quantitative  SIMS  data.  From  calculations,  and  by  comparison  to 
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Figure  23  -  Dependence  of  the  Hl+,  Sly-*- ,  and  Xl3+  Ion  yield  (most 

Intense  Isotope  peaks  at  m/e  *  58,  116,  174)  on  Oy  exposure. 

Sly"1"  and  Nlj'*"  Intensities  are  corrected  for  Isotope 
distribution  and  quadrupole  transmission.  The  corrected 
Sl3+  Intensity  Is  magnified  by  a  factor  of  3.  From  reference 
106. 

similar  studies  by  Oechsner  on  oxidized  Ce  and  Ta  (107,108),  we  anticipate  that  the  Ml 
neutral  yield  would  decrease  by  at  most  a  factor  of  2  when  the  oxygen  Is  adsorbed.  Similar  - 
observations  for  negative  Ion  enhancements  have  been  made  by  Yu  for  Cs  adsorption  on 
various  metals  (109).  Thus,  oxygen  which  often  Increases  the  the  surface  work  function 
enhances  positive  Ions,  while  Cs  which  usually  decreases  the  work  function,  enhances 
negative  Ions,  Neither  of  these  atoms  has  a  major  effect  on  neutral  yields.  This 
observation  has  been  used  for  many  years  with  the  Ion  microprobe  since  sensitivities 
(110,111)  could  be  greatly  Increased  by  bombarding  with  0y+  or  Cs+,  presumably  since  the 
primary  Ions  subsequently  reacted  with  atoms  In  the  near  surface  region  to  form  new 
molecular  species.  So  as  not  to  present  a  completely  consistent  picture,  however.  It  should 
also  be  pointed  out  that  Oy4,  bombardment  has  been  observed  to  enhance  negative  ions  (112). 

For  a  given  single  crystal  matrix,  the  Ionization  probability  has  been  found  to  be 
reasonably  Independent  of  the  azimuthal  ejection  angle.  For  Cu  single  crystals  Yurasova  and 
co-uorkers  (113)  found  intensity  variations  with  crystal  rotation  nearly  Identical  to  those 
reported  by  Wehner  for  Cu  neutrals  (15).  Similar  observations  have  been  reported  for  Si 
(114).  The  polar  angle  distributions  for  Sl+  are  also  quite  similar  to  those  calculated  for 
N1  atoms  If  an  image  force  of  3.6  eF  Is  allowed  to  bend  the  neutral  trajectories  away  from 
tha  surface  normal  <115  ).  These  effects  will  be  discussed  In  more  decail  later.  We  should 
also  note  that  rather  strong  variations  In  neutralization  probability  with  angle  have  been 
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noted  when  the  primary  Ion  directly  scatters  from  the  surface.  For  example,  by  comparing 
che  scattered  He*  intensity  to  the  scattered  Li*  intensity  at  600  eV,  Taglauer  and 
co-workers  found  chat  azimuthal  rotation  of  a  Ml( 1 10)  crystal  did  not  alter  che  He* 
intensity  (116).  The  Li*  beam,  which  is  never  neutralized,  exhibited  large  fluctuations  due 
to  second  layer  scattering.  Thus  they  concluded  that  when  the  Re*  penetrated  below  the 
first  layer,  it  was  effectively  neutralized  but  during  surface  collisions,  neighboring  atoms 
did  noc  Influence  che  ion  yield.  Woodruff  and  co-workers  have  scudied  angular  variations  of 
He*  from  Si(001)  covered  with  25Z  of  a  monolayer  of  oxygen  (p(2x2)  overlayer)  (117,118). 

They  found  chat  when  the  He*  beam  scattered  near  a  protruding  surface  atom  adjacent  to  che 
Impact  point  chat  it  experienced  a  higher  neutralization  probability.  At  this  point,  then, 
there  is  conflicting  evidence  as  to  how  similar  atom  yields  and  ion  yields  are  as  a  function 
of  angle. 

The  last  observable  property  of  importance  is  the  secondary  particle  energy 
distribution.  Qualitatively,  che  ions  and  neutrals  have  similar  energy  distributions  in 
that  ions  with  energy  ranging  from  nearly  zero  to  several  hundred  eV  have  been  observed. 

More  specifically,  che  ion  distribution  tends  to  be  somewhat  broader  than  the  neutral 
distribution  and  tails  off  as  E”13*5  rather  than  as  E~z -Q  as  for  neutrals.  A  detailed 
comparison  has  been  made  between  Cu*  ions  and  Cu  atoms  ejected  from  a  polycrystalline  Cu 
surface  (2d)  as  shown  in  Figure  2d.  Similar  trends  have  been  observed  for  many  systems.  One 
known  exception  to  this  observation  is  that  low  ionization  potential  elesiencs  like  Na*  and 
K*  appear  to  have  energy  distributions  either  Identical  to,  or  narrower  than,  che 
corresponding  neutrals.  This  has  been  observed  for  He*  in  silica  (119)  and  for  K*  on  Hi 
(120).  In  any  case,  however,  there  does  not  appear  to  be  a  strong  dependence  of  R*  on 
particle  velocity  since  the  energy  distributions  of  ions  and  neutrals  are  at  least 
qualitatively  similar. 

As  indicated  above,  there  have  been  many  attempts  to  explain  some  or  all  of  the  above 
differences  between  ions  and  neutrals.  In  this  review,  all  of  these  models  will  not  be 
recalled,  racher  we  will  focus  on  some  Important  historical  developments  and  comment  on  some 
of  the  moat  promising  new  approaches  to  straightening  out  che  mess.  Of  particular  lncerest 
will  be  chose  models  chat  take  into  account  the  complexity  of  che  sputtarlng  evenc  and 
actempt  to  understand  ionization  from  a  microscopic  point  of  view. 

Probably  the  most  significant  initial  observation  regarding  che  ionization  probability 
was  that  the  value  of  R*  depends  exponentially  on  ionization  potential  t  as 

R+  a  sxp(-I/K)  (35) 

where  X  is  a  fitting  parameter  (105,121).  The  approximate  nature  of  this  relationship  is 
illustrated  in  Figure  25  for  a  nunber  of  the  elements  (121).  Although  che  form  of  Souation 
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.2*  "  {*)  TfPlcal  energy  spectre  of  sputtered  Cu  neutrals,  obtained 
by  electron  bombardment  Ionization. 

<b)  Typical  energy  spectra  of  sputtered  Cu  lone,  showing 
a  broader  peak  than  that  of  the  neutrals  tn  (a).  From 
reference  24 . 

35  seems  as  good  as  any  considering  the  experimental  difficulties  In  measuring  ,  it  has 
been  a  major  chore  to  relate  the  result  to  a  physically  reasonable  model.  In  a  widely 
criticized  report.  Andersen  and  Hlnthorn.  U21)  utilized  the  Saha  equation,  which  are 

derived  from  Boltzmann  statistic,  applied  to  a  plasm,  equilibrium  of  element  M  with  its 
ionized  species  as 


+  It' 


(36) 


SPUTTEREO  POSITIVE  IONS 


CUNOPYROXENE  (AC-362-C1) 
“O",  16hV 


(a)  (b) 


Figure  25  -  (a)  Saha-Sggert  loolzadon  equation  applied  to  sputtered 

nsltlve  ion  data  from  a  mineral  standard  bombarded  with 
0-  at  16  IcV. 

(b)  Saha  equation  applied  to  sputtered  negative  Ion  data 
from  a  mineral  standard  bombarded  with  160“  at  11  '*V.  From 
reference  121. 


The  resulting  equilibrium  expression  is  then 

Cjj+CV/Cm’  -  [(Zma+mekD^/Z/hmol-Ca+JSe/BoJexpC-I/kT)  (37) 

where  the  3's  are  Che  appropriate  partition  functions,  and  the  C's  represent  che 
concentration  of  che  Indicated  species.  The  temperature,  T,  Is  essentially  a  fitting 
parameter  for  a  series  of  elements  ejected  from  a  common  matrix.  Typical  values  of  T  range 
from  2000  to  15000  K.  There  are  many  reasons  why  this  explanation  of  Equation  35  is 
unphyslcal  but  perhaps  che  most  persuasive  is  from  results  of  the  molecular  dynamics 
calculations.  These  calculations  show  chat  equilibrium  Is  not  attained  anywhere  near  che 
target  atom  during  particle  ejection,  a  result  graphically  depicted  in  Figure  1.  Further, 
the  energy  distribution  of  the  ejected  particles  should  have  a  Maxwell-Boltzmann  shape,  a 
result  not  predicted  by  collision  cascade  models  or  by  axperlmental  results. 

A  better  explanation  of  Equation  35  as  proposed  by  Williams  (122,123)  is  that  che 


sputtered  atom  Interacts  electronically  tilth  the  atoms  near  the  ejection  site  via  resonant 
electron  exchange  processes.  Thus,  if  an  energy  level  in  the  surface  is  occupied  by  an 
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electron  and  an  equivalent  energy  level  Is  unoccupied  in  the  ejecting  atom,  electron 
exchange  nay  occur.  This  situation  can  be  approximated  by  Equation  35  assuming  that  the 
excited  electrons  are  at  equilibrium  and  nay  be  described  using  Ferml-Dirac  statistics. 
Specifically,  the  model  assumes  that  the  yield  of  positive  ions,  Y*  is 

T+  •  TP,  (38) 

where  Pe  is  the  probability  that  the  electronic  level  in  the  surface  at  energy  1  is  vacant. 
The  value  of  Pe  may  be  defined  as 

Pa  -  U+exp[(I-Et)/kTe] ("I  (39) 

where  £c  is  the  highest  occupied  state  in  the  quiescent  surface  and  Ta  is  the  equilibrium 
temperature  of  the  electrons.  In  a  metal  Et  corresponds  roughly  to  the  work  function  of  the 
surface  and  in  an  lnsulacor  or  semiconductor  to  the  valence  band  edge.  Since  the  yield  of 
atoms  ejecting  in  a  neutral  state  is 


Y*  -  YPf 

(40) 

and 

Pe  +  Pf  -  1 

(41) 

Chen 

g+  .  y«7Y*  -  exp[-(  I-Et)/kT,] 

(42) 

The  schematic  electronic  density  of  states  and  the  potential  energy  diagram  are  Illustrated 
in  Figure  26. 


Surface  Atom 


Electron  Energy  Potential  Energy  Diagram 

Distribution 


Figure  26  -  Surface  electron  excitation  model.  From  reference  122. 
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Although  this  model  remove*  many  of  tha  objections  noted  for  the  Andersen-Hlnthorne 
rationalization.  It  Is  not  readily  expanded  to  Include  an  explicit  representation  of  the 
bulk  (or  surface)  density  of  electronic  states  nor  to  Include  the  dependence  of  R+  on  angle 
or  energy.  (There  Is  an  iapllclc  approximate  exp(-conat/»2)  dependence  of  R+  on  ion 
velocity  In  Equation  42  as  discussed  in  Reference  122).  To  incorporate  these  aspects  of  the 
problem  It  Is  necessary  to  know  something  about  the  wavefunctlons  of  the  atoms  In  the  solid 
and  co  relate  these  wavefunctlons  to  electron  exchange  probabilities  between  two  Interacting 
systems. 

Sroubek  and  co-workers  have  made  considerable  progress  along  these  lines  (124).  Their 
first  approach  Is  to  consider  a  single  time  independent  level  of  an  adsorbate  atom  $a  of 
energy  ea  that  interacts  with  a  single  time  independent  metal  substrate  level  of  energy 
£i .  The  electron  is  assumed  co  reside  initially  in  ta  on  the  adatom  with  level 
unoccupied.  If  the  time  dependent  coupling  (hopping  Integral)  between  Che  levels  V(t)  Is 
slowly  turned  on,  then  the  R4-  can  be  calculated  from  perturbation  theory  as  (125,126) 

H4"  -  |  l/h7v(t)expfl(ei-sa)t/* ]dt  |  2  (43) 

Many  functional  dependences  for  V(t)  have  been  tested  although 

7(0  -  V0*xp(-lvt),  t>0  (44) 

gives  the  most  physically  realistic  values  of  R*.  Here,  1  Is  a  constant  which 
describes  the  overlap  of  wavefunctlons  and  Is  approximately  2  .  It  is  of 

interest  chat  If  the  two  atoms  move  apart  with  constant  velocity  v  from  a 
starting  separation  ra,  the  coupling  can  be  rewritten  as 

V(t)  -  V0*xp{-Xtr(t)-r«H  (45) 

Solving  Equation  43  using  V(t)  from  Equation  44  yields  for  the  ionization 
probability 

R4-  -  «2V02v2X2/(e1-ea)A  (46) 

Mote  that  H4-  depends  on  the  square  of  the  particle  velocity.  The  modal  also 
predicts  strong  dependence  on  the  coupling  between  Che  adatom  and  Che  surface 
acom,  both  the  strength  V„  and  range  1,  and  even  stronger  sensitivity  to  the 
relative  position  of  the  electronic  levels. 

Although  this  approach  accounts  for  tha  detailed  electronic  structure  of  the  relevant 
atoms.  It  clearly  Is  unrealistic  because  only  two  levels  are  Involved  In  the  collision.  It 
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Is  possible,  however,  Co  generalize  Equation  43  such  chat  it  Includes  a  continuum  of  levels 
and  allow  for  the  tunneling  of  electrons  froa  occupied  scaces  below  the  Feral  level  or  froa 
Auger  processes  027).  This  expression  for  K+  nay  be  obtained  as 

-  2/ir{exp[-C1»(I-e)/hXvlt  (47) 

using  che  coupling  Integral  given  In  Equation  44.  Here,  Cj  Is  a  constant  and  4  Is  the  work 
function  of  che  surface.  Note  chat  this  result  very  closely  resembles  che  original 
Andersen-Hlnchorne  observation.  The  result  differs  quite  draaatlcally  from  Equation  46  in 
that  It  exhibits  no  dependence  on  the  coupling  strength  7„  and  has  a  very  strong  dependence 
on  velocity  (exp(-consc/v) ) .  A  major  difficulty,  however.  Is  that  the  effective  temperature 

«  Xvh/xk  (48) 

Is  much  smaller  than  experimentally  observed  (127,128)  and  calculated  values  of  HT*"  are  many 
orders  of  magnitude  too  small. 

It  Is  of  Interest  that  a  still  different  velocity  dependence  may  be  predicted  by 
assuming  that  the  adatom  level  ea  changes  In  dme  (or  distance)  as  che  particle  leaves  the 
Influence  of  the  band  of  levels  In  the  metal  (125).  This  corrected  expression  Is 

a  [(ea-Ef)Xv’l/2/AE  (49) 

where  Ef  Is  the  feral  energy  and  AE  Is  the  width  of  the  conduction  band.  Thus,  there  are  a 
range  of  models  which  predict  the  dependence  of  R*  on  velocity  to  vary  from  v»/2  to  v^  to 
exp(-const/v)  to  exp(-const/v2) . 

In  boch  the  two  level  model  and  the  model  where  a  continuum  of  energy  levels  are 
Included,  che  explicit  electronic  structure  of  che  ejecting  atom  and  the  solid  are  still  not 
considered  In  any  detail.  Another  difficulty  -  and  one  which  Is  Inherent  In  virtually  every 
model  discussed  so  far  -  Is  chat  che  electronic  structure  of  the  solid  must  surely  change  as 
che  collision  cascade  evolves  In  time.  For  example,  the  electronic  structure  of  a  metal 
crystal  consisting  of  a  periodic  array  of  atoms  can  be  described  quantum  mechanically  using 
che  appropriate  Bloch  wavefunctlons  and  eigenstates.  After  che  Impact  of  the  primary  ion, 
however,  this  periodicity  la  removed  completely  as  the  atoms  more  closely  resemble  a  liquid. 
In  fact  aa  has  been  pointed  out  by  Williams  (122),  since  the  atoms  are  moving  with  energies 
of  several  eV  per  atom,  acom-atom  distances  can  vary  over  a  much  wider  range  chan  normally 
found  In  amorphous  materials.  This  disorder  may  have  che  effect  of  localizing  the  Itinerant 
electrons  of  a  metal  or  of  creating  new  states  with  energies  within  che  bandgap  of 
semiconductors  or  Insulators.  This  effect  is  illustrated  in  Figure  27.  Thus,  it  may  turn 
out  chat  any  model  which  relies  on  che  band  structure  of  the  unperturbed  solid  to  calculate 


semiconductors  as  a  function  of  energy,  (a)  crystalline 
semiconductor,  (b)  amorphous  but  stoichiometric  semiconductor, 

(c)  probable  density  of  state  for  highly  disordered 
semiconductor.  From  reference  122. 

R+  already  has  Inherently  gross  errors.  What  Is  really  required  is  a  formalism  where 
electronic  motion  and  nuclear  motion  are  determined  in  time  together.  This  approach  has 
been  attempted  for  a  few  simple  situations  and  offers  considerable  promise,  not  necessarily 
to  generate  quantitative  predictions  of  RT*" ,  but  to  provide  insight  into  the  important 
factors  which  Influence  surface  ionization. 

To  consider  both  electronic  and  nuclear  motion,  we  begin  using  a  Hamiltonian  which  is 
separated  into  a  part  chat  describes  the  motion  of  the  atomic  nuclei  and  a  part  that 
describes  the  motion  of  the  electrons.  The  nuclear  positions  are  determined  using  the 
procedure  outlined  earlier  in  this  seccion  and  provide  the  appropriate  tine  dependent 
coordinates  for  solution  of  the  electronic  part  of  the  problem.  This  approach  presumes  that 
there  is  no  coupling  between  electronic  and  nuclear  motion  and  that  Che  classical  trajectory 
is  unaffected  by  the  electronic  transitions  that  occu  during  the  trajectory.  The 
conditions  where  this  assumption  is  valid  have  been  outlined  In  some  detail  by  Tully  (129). 

To  determine  H*  for  an  adatom  ajecte.  from  a  surface,  it  is  necessary  to  solve  the 
time  dependent  Schrodinger  equation  for  the  system  using  the  time  dependent  electronic 


Hamiltonian 


Here,  Hst  Is  an  atomic  Hamiltonian,  and  0  is  i  tima  or  position  dependent  inter-atomic 
electronic  interaction  integral.  The  coupling  interaction  neglects  electron-electron 
interactions  so  no  two  electron  or  Auger  type  processes  are  considered.  The  solutlona  to 
Hat  are  given  by 


9it*e  "  (51) 

where  the  4y' i  are  the  orthonormal,  atomic  vavafunctions  for  Che  atoms  in  the  solid. 

To  solve  the  Schrodlnger  equation  using  the  Hamiltonian  in  Equation  SO  we  then  seek 
time  dependent  molecular  vavafunctions  where 

)*i 

Hilt  -  tfi  —  (52) 

Although  there  are  number  of  possible  approaches  to  representing  the  efforts  so 
far  have  focussed  on  utilizing  the  tight  binding  approximation.  In  this  case,  the 

electronic  levels  in  the  solid  may  be  expresssed  as  a  linear  combination  of  atomic  orbitals 
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’l  ’  ^cvi(t)exp(-iElt/4>)  (53) 

where  an  arbitrary  phase  factor  has  been  Included.  The  expansion  coefficients  are  time 
dependent  since  the  position  of  the  nuclei  will  be  changing  during  the  evolution  of  che 
collision  cascade.  Combining  Equations  52  and  53  integrating  over  the  electronic 
coordinates  yields  the  following  sec  of  coupled  equations: 

ifjdcki.(t)/dt  1  *  U  Ui>cij(t)  +  (  efc-Sj^Cfci^t)  (56) 

These  equations  of  motion  for  the  Cfci( O ' s  can  be  Integrated  simultaneously  with  the 
classical  equations  for  the  nuclear  motion.  As  was  used  in  che  evo  level  case,  the  coupling 
matrix  element  or  hopping  integral  is  given  by 

<«t|u|$3>  »  v[r(t)1  •  Vfcj  exp[“lkj(rkj-r*Vtj5 1  ;*lcj  4  0 ,ic*j  (55) 

-  0  jritj  -  O.Mj  (56) 

The  value  of  may  be  calculated  from  che  coupling  coefficients  in  che  following 
manner : 

1)  Choose  soma  atomic  energies  and  parameters  for  7[r(t)|. 

2)  Solve  for  the  molecular  energies  E^  at  time  t»0  by  setting  Equation  56  equal  co 


zero. 


3)  Determine  the  nuclear  motion  of  Che  total  systea  by  the  molecular  dynaalcs  method. 

At  each  tloestep,  Che  new  nuclear  poaitlona  serve  as  Input  to  Che  coupling  matrix 
eleaent.  The  coupled  equations  of  aoclon  for  the  cvi(c)'s  are  simultaneously 
integrated  with  the  nuclear  poaitlona  and  velocities. 

4)  Evaluate  at  the  end  of  cha  collision  process  by  projecting  the  adatoa  state 

on  all  states  chat  were  originally  unoccupied.  FT  is  thus  given  by 


unocc 

**■  -  I  |cai(t-«)J2  (57) 

The  numerical  analysis  for  this  prescription  is  currently  beyond  the  rela  of  aosc 
present  day  computers.  Currently,  Che  solid  has  been  approxiaated  using  a  linear  chain  of 
atoms  (130)  and  by  a  cluster  of  five  atoms  with  the  adsorbate  atom  placed  on  the  surface 
(129)  as  shown  in  Figure  28.  Another  approach  involves  monitoring  the  nuclear  aoclon  of  a 
large  collection  of  particles  (~300  or  sore)  while  only  selecting  a  few  (5-10)  to  folio* 
electronically  (131). 


Figure  28  -  The  atomic  model  used  in  the  computer  simulation.  The  model 

consists  of  the  atomic  substrate  (white  spheres),  the  particle 
which  is  sputtered  (shaded  sphere),  and  the  bombarding  parcicle 
(black  sphere).  The  Interatomic  distances  are  marked  in  4.  The 
electronic  structure  of  the  model  is  shown  in  the  inset  with 
electrons  depicted  by  black  dots.  1  is  the  ionization  energy 
of  the  sputtered  particle  and  is  assumed  co  be  S.25  eV  in  this 
case,  t  marked  the  highest  density  of  the  empty  electronic 
states.  Froa  reference  128. 


The  incorporation  of  the  nuclear  motion  into  the  ionization  theory  is  critical  co 


obtain  physically  reasonable  results.  In  Sroubek's  model  study,  he  found  chat  by  varying 
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che  ionization  energy  of  che  ejecting  adsorbate  atom,  he  could  produce  values  of  R*  which 
could  be  fie  to  Equation  47  with  effective  temperatures  between  1800  and  3300  K.  He  could 
also  obcain  EC*  values  that  ranged  from  10” *■  to  10~$,  much  more  reasonable  values  chan  are 
obealned  from  lattices  thee  are  not  allowed  to  move  (128,130).  Also,  of  interest  is  that 
for  his  initial  studies,  EC*  exhibited  a  very  weak  dependence  on  particle  velocity.  Although 
che  development  of  this  procedure  is  in  its  very  esrly  stages,  it  appears  that  all  of  che 
fundamental  ingredients  are  available  to  construct  a  general  understanding  of  che  ionization 
phenomenon. 

There  are  many  other  theories  chat  have  been  sec  forth  to  try  to  predict  R*,  although 
virtually  all  of  these  are  qualitative  In  nature.  For  example,  in  the  'bond-breaking-  model 
(132)  It  is  proposed  chat  if  a  molecule  is  already  charged  in  the  original  solid  chat  the 
atomic  component  may  be  ejected  with  its  charge  state  intact.  The  high  yield  of  Ha*  and  Cl~ 
from  HaCl  illustrates  this  point  nicely.  Other  models  that  Incorporate  surface  work 
functions  (133)  or  bulk  band  structure  (134)  to  predict  EC*  suffer  from  the  fact  that  these 
parameters  have  little  meaning  during  che  collision  cascade. 

Anocher  important  theoretical  approach  focusses  on  the  fact  that  once  the  ion  is 
formed  during  che  excitation  process,  it  will  undergo  de-excltatlon  as  long  as  it  is  close 
to  che  surface  (135).  Two  types  of  neutralization  have  been  proposed  for  ions  approaching  a 
metal  surface.  The  first  type  involves  resonance  tunneling  of  an  electron  from  che  filled 
levels  of  che  metal  to  an  excited  level  of  Che  ion.  The  second  type  Involves  two  electron 
processes  where  an  electron  from  the  metal  fills  a  lower  energy  level  of  the  ion.  The 
excess  energy  is  lost  via  emission  of  an  Auger  electron  from  either  the  metal  or  the  ion. 

The  transitions  are  illustrated  schematically  in  Figure  29.  The  probability  of  a  particle 
escaping  de-excitation  is  given  as 


P( s-“)  -  axp£-A/av)  (58) 

where  s  is  the  distance  from  the  surface,  A  is  che  transition  race  of  an  electron  from  the 
surface  to  che  ion,  i  is  a  critical  distance  and  v  is  che  particle  velocity  normal  to  the 
surface . 

The  energy  distribution  of  sputtered  ions  should  have  a  similar  velocity  dependence  as 

y*(v)  »  y*(v)exp(-A/av  )  (59) 

where  y*(v)  may  be  estimated  from  Equation  7  or  ocher  models  available  for  predicting  the 
velocity  dependence  of  the  neutrals.  Whan  Equation  59  is  fit  to  che  energy  distribution  of 
Cu  atoms  (136)  a  value  of  A/a  of  2.5  x  10^  to  1.7  x  10*  cm/ s  is  obtained  which  is  comparable 
to  values  of  A/a  used  to  explain  photon  emission  data  (137). 
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Although  this  development  seems  physically  reasonable  from  a  number  of  viewpoints , 
there  Is  a  major  difficulty  In  applying  Equation  59  to  any  real  situations.  The  treatment 
assumes  that  the  lattice  effectively  remains  static  during  sputtering.  As  discussed 
earlier,  Sroubek  has  shown  thac  the  lattice  motion  may  effect  both  the  total  yield  y-t  and 
the  velocity  dependence  of  T*-  in  a  major  fashion.  In  fact.  It  appears  chat  by  using  the 
coupling  Integral  given  In  Equation  14  that  Sroubek's  approach  and  the  Ion  neutralization 
approach  will  yield  similar  results  In  the  limit  of  a  static  lattice.  Thus,  Sroubek  has 
noted  thac  Equation  59  predicts  a  velocity  dependence  thac  Is  much  coo  strong  (128). 

There  are  very  few  experimental  studies  which  have  been  capable  of  helping  to  decide 
which  theories  of  Che  Ionization  process  ara  Indeed  the  most  reliable.  From  a  quantitative 
point  of  view,  however,  the  velocity  dependence  of  K*  varies  considerably  from  model  to 
model  and  would  therefore  be  a  valuable  functionality  to  obtain  experimentally.  There  have 
been  two  quite  detailed  scudles  chat  have  recently  been  completed  chat  qualitatively  support 
Sroubak's  contention  chat  the  velocity  dependence  of  R*-  should  not  be  particularly  strong. 

first  study  focusses  on  Che  adsorption  of  CO  on  Sl(OOl)  since  this  system  produces 
Intense  Nl*  signals  and  since  the  surface  adsorbate  structure  has  been  well  characterized  by 


ocher  techniques  (113).  The  measurements  were  performed  using  a  new  SIMS  apparatus  aimed  at 
determining  ion  yields  at  all  9  and  »  ejection  angles  for  a  normally  Incident  beam  (120). 

As  shown  in  Figure  30  chis  measurement  was  carried  out  by  utilizing  a  quadrupole  mass 
spectrometer  chat  could  be  rotated  under  tJHV  conditions  with  respect  to  Che  incident  ion 
beam.  In  addition,  a  medium  resolution  energy  selector  consisting  of  a  90*  spherical  sector 
provided  the  capability  of  examining  how  the  angular  distributions  changed  with  secondary 
ion  energy.  The  idea  behind  designing  chis  apparatus  was  that  the  ion  trajectories  could  be 
directly  compared  co  those  calculated  for  the  neutrals  using  the  molecular  dynamics 
treatment . 


Include  M,  crystal  manipulator;  Q.M.S.,  quadrupole  mass 
spectrometer;  I.G.,  primary  ion  source;  E.S.,  energy 
spectrometer;  G,  Bayard-Alpert  gauge;  T,  crystal  target;  and 
G.I.,  gas  Inlet.  Auxiliary  components  are  omitted  for 
graphical  clarity.  The  SIMS  experimental  geometry  and  coordinate 
system  are  defined  in  the  inset.  From  reference  120. 

The  results  of  the  measurements  have  so  far  only  been  analyzed  for  the  MI  species  due 
to  the  difficulty  in  obtaining  enough  particles  in  the  calculation  that  are  ejected  at 


specific  angles  and  energies.  There  have  been,  however,  several  Important  consequences  of 
the  experiments.  First,  If  the  CO  Is  assumed  to  bond  In  the  atop  site,  the  polar  angle 
distributions  of  Ml'*'  thac  are  measured  experimentally  do  not  agree  well  with  those 
calculated  for  a  Nl(OOl)  surface  covered  with  501  of  a  monolayer  of  CO.  As  shown  In  Flgui t 
31,  for  example,  calculated  polar  plots  peak  at  an  angle  of  about  10*  closer  to  the  normal 
than  the  measured  ones  (11S).  Thus,  If  It  Is  assumed  that  the  calculations  are  correct, 
there  must  he  a  fundamental  difference  between  Ion  trajectories  and  neutral  crajectorlee. 
This  disci epency  has  been  resolved  by  Incorporating  the  presence  of  the  classical  Image 
force  Into  the  computational  model.  If  It  is  assumed  chat  the  particle  Instantaneously 
becomes  an  Ion,  E<„.g,.  Is  given  by 

Elmage  “  e2/4*  *  3.6  eV/a0  (60) 

where  ac  Is  the  height  In  1  of  the  particle  above  the  jelllum  step-edge  at  the  Instant  of 
Ionization.  If  9n  Is  Che  angle  of  the  atom's  velocity  vector  at  this  point,  the  ion  emerges 
with  a  final  direction  given  by 

*1  -  tan"1[(E0sin29n)/(H0coa29n-EiBagB>  )1/2  (61) 

where  9^  Is  che  corrected  polar  ejection  angle  as  measured  from  the  surface  normal  of  the 
leaving  Ion,  and  E0  Is  the  kinetic  energy  of  the  neutral  particle  ejected  ac  a  polar  angle 
®n*  The  corrected  theoretical  energy  distributions  using  Elmage  *  3 .6  eV  are  also  shown  In 
Tlgure  31.  The  agreement  with  experiment  over  all  polar  angles  Is  now  quite  good.  Note 
also  that  che  width  of  the  polar  angle  distribution  is  much  sharper  ac  higher  kinetic 
energies.  This  expected  result  occurs  since  the  particles  ejected  early  In  che  collision 
cascade  see  more  surface  order  than  chose  ejected  lace,  after  the  crystal  structure  Is 
nearly  destroyed. 

To  check  chat  che  success  of  this  simple  correction  Is  not  just  coincidence.  It  Is 
possible  to  compare  both  che  azimuthal  distributions  and  che  energy  distributions  of  the 
measured  Nl+  Intensity  to  calculated  values.  For  each  comparison,  the  agreement  has  been 
very  good  using  the  same  value  of  Elmage  (115).  In  fact.  It  haa  been  proposed  that  the 
reason  chat  che  energy  distribution  for  ions  Is  usually  broader  chan  for  neutrals  Is  due 
solely  to  the  presence  of  Eim^gg  (138).  In  relation  to  Thompson's  energy  distributions,  the 
Image  potential  acts  as  an  additional  surface  binding  force  such  chat  Equation  7  becomes 


Y*(E)aCE/(E+Eb-EiMg,)2 


(62) 
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Figure  31  -  Polar  angle  distributions  for  N1  ejected  from  Si(001)c(2x2)-C0 
The  upper  set  of  curves  Is  recorded  for  a  secondary  Ion  energy 
of  712  eV  while  the  lower  sec  of  curves  is  taken  at  2212  eg. 

From  reference  U5. 

A  further  point  Is  chat  the  classical  dynamics  model  assumes  chat  there  Is  no  velocity 
dependence  or  azimuthal  dependence  to  if .  Since  this  model  can  accurately  describe  the 
experimental  situation.  It  must  be  concluded  chat  In  this  case  these  properties  of  R+  are  at 
least  approximately  valid. 

In  a  completely  different  approach,  Tu  (139)  has  Investigated  the  velocity  dependence 
of  R+  by  assuming  that  the  functional  form  of  the  Jforskov  and  Lundquist  prediction  Is  Indeed 
correct.  By  rewriting  Equation  47  as 

R+  -  (2/v)expr(p-I)/e0;  (63) 

where  eg  •  Ci*,  Tu  reasoned  that  ac  constant  v^,  a  plot  of  log  fT*"  va  (g-I)  should  yield 

a  straight  line  with  slope  l/e0.  At  different  v^,  however,  the  slope  should  vary  according 
to  the  definition  of  e0.  He  was  able  to  vary  )  by  a  sufficient  amount  by  evaporating  very 
small  quantities  of  LI  onto  the  sample  surface.  The  results  for  0"  emission  from  Ion 
bombarded  oxidized  vanadium  are  shown  In  Figure  32.  The  slope  of  each  line  should  give  the 
value  of  e,. _ When  the  e0  values  are  then  plotted  vs  vcoad  as  shown  In  Figure  33,  the  values 


59*  EMISSION 


Figure  32 


Figure  33 


The  work  function  dependence  of  Che  0~  yield  at  two  different 
emission  energies  8.3  and  63  eV  but  for  the  same  angle  of 
emission  35* .  Notice  the  exponential  dependence  on  do  for  -do 
<2  eV.  From  reference  139. 


The  dependence  of  e„  for  0”  sputtered  fron  oxygenated  (1U 
vanadium  surfaces  on  the  normal  component  v  of  the  emission 
velocity.  The  approximate  proportionality  between  e0  and  v 
for  Vj.  >1x10®  cm/ sec  Is  In  agreement  with  theory.  From 
reference  139. 
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should  fall  on  a  scralghc  llna  as  predicted  by  Equation  47.  Tha  results  show  that  at 
■clnetlc  energies  less  than  about  20  eV  that  e0  shows  alisost  no  velocity  dependence  but  that 
[  at  energies  greater  than  20  eV,  Equation  47  becones  valid.  Since  the  peak  in  the  energy 

distribution  occurs  at  leas  than  10  eV ,  these  experlnenta  clearly  indicate  that  T*"  also 
exhibits  a  weak  velocity  dependence. 

We  have  concentrated  in  this  discussion  on  mechanisms  of  ionization  that  occur  in 
concert  with  the  development  of  che  collision  cascade.  Other  possible  processes  are  thought 
to  exist  which  lead  to  ion  ejection.  For  example,  Williams  has  proposed  that  since  positive 
halogen  ions  eject  with  an  energy  distribution  which  is  only  abouc  15  eV  wide  and  peaks  at 
6-8  eV,  che  ion  is  noc  produced  by  direct  colllsional  events  (140).  The  ionization  is 
proposed  co  proceed  via  Interatomic  Auger  decay  followed  by  coulomblc  ejection  of  the  ion  in 
a  mechanism  similar  co  thac  proposed  by  Knotek  and  Felbelman  to  explain  electron  and  photon 
stimulated  desorption  (141).  dearly,  none  of  the  above  models  would  be  appropriate  for 
this  situation. 

To  conclude  this  section,  then,  it  is  clear  chat  much  work  as  been  devoted  co 
understanding  the  ejection  of  neutrals  and  ions  from  ion  bombarded  surfaces.  Our 
understanding  of  the  neutral  ejection  process  is  somewhat  more  advanced  chan  for  che  Ions. 
However,  there  is  now  a  formalism  available  for  coupling  che  microscopic  nuclear  moclon 
surrounding  che  ejection  site  co  che  possible  electron  transfer  processes  thac  can  occur  in 
the  solid.  For  all  of  che  cheorles,  there  are  many  predictions  chat  can  be  experimentally 
tested.  In  che  next  sections  we  will  attempt  co  build  upon  the  ideas  developed  by  the 
theorists  to  attempt  to  mors  fully  understand  the  resulting  mass  spectra  and  to  see  if  we 
can  learn  something  abouc  che  structure  and  chemistry  of  the  sample  itself. 

IV.  Structure-Sensitive  Factors 

In  previous  sections,  there  .ias  been  a  qualitative  description  of  a  varlecy  of 
spuccerlng  experiments  that  clearly  contain  information  about  the  crystallography  of  the 
solid.  Uehner's  spots  are  probably  the  mosc  spectacular  examples  of  this  fact.  In 
addition,  it  appears  thac  the  state  of  the  theory  which  describes  che  Ion  impact  event  is 
sufficiently  advanced  such  that  a  oora  quantltatlva  understanding  of  che  scructure-senslclv 
factors  In  sputtering  should  be  forthcoming.  Through  these  comparisons,  it  may  also  be 
possible  co  use  Ion  beam  methods  to  learn  something  new  about  the  surface  structure  of 
previously  uncharacterized  materials.  In  this  section,  then,  mt  will  examine  these 
structure  sensitive  factors  for  che  ejecclon  of  atoms  (or  Iona)  from  single  crystal  samples 
and  sea  how  well  the  measured  observables  agree  with  che  choeen  theory. 


A.  Yields 


Recall  first  chat  the  total  yield  of  pertlclea  ejected  froa  a  copper  single  crystal 
surface  Is  a  function  of  cryscal  orientation  with  Y(  1 ll)>Y(001)>Y( 110) .  The  experiments  to 
determine  this  relation  (18)  were  performed  uelng  the  weight-loss  technique  and  are 
therefore  subject  to  a  certain  degree  of  uncertainty  becauee  of  surface  damage  by  the 
primary  ion.  In  addition,  Che  vacuum  conditions  were  rather  poor  by  today's  standards.  Yet, 
chese  observations  have  never  been  contradicted,  although  It  would  be  extremely  helpful  to 
repeat  many  of  them  using  as  low  a  primary  dose  as  possible  and  using  simultaneous  LEED 
monitoring.  Ic  Is  possible  to  examine  the  reasons  for  chese  yield  differences  using  the 
molecular  dynamics  calculations.  These  calculations  were  first  performed  only  recently  on 
Cu  since  large  model  crystals  are  needed  In  order  to  sufficiently  contain  the  nuclear 
motion  chat  leads  to  particle  ejection.  The  required  computer  time  Is  nacher  lengthy,  even 
for  the  fastest  computer  systems.  The  results  of  these  calculations  are  given  In  Table  II 
(142).  Note  thac  Che  relative  order  of  the  calculated  yields  agree  with  that  obtained  from 
the  early  experiments.  A  detailed  analysis  of  the  Individual  trajectories  provides  some 
Insight  into  why  this  order  is  observed.  In  general,  Che  momentum  resulting  from  the 


Table  II 

Sputtering  Yield  of  600  eV  Argon  on  Copper( 100) ,(110)  and  (111) 


Cryscal 

face 

Total  Number 
of 

trajectories 

Total  Number 
of  atoms 
ejected 

Sputtering 

Yield 

Relative  Yield 

Calculated8  Experimental b 

(100) 

111 

436 

3.93 

1.00 

1.0 

(110) 

121 

429 

3.54 

0.90 

0.5 

(111) 

108 

700 

6.48 

1.65 

1.3 

^Normalized  to  the  (100)  face. 

^Measured  at  1  <eV  by  weight  loss,  normalized  to  the  (100)  face  (18). 

primary  Ion  Is  kepc  close  to  the  surface  If  the  density  of  atoms  Is  high  as  In  the  case  for 
Che  (111)  orientation.  For  much  more  open  configurations  such  as  che  (110)  orientation  Che 
momentum  more  readily  penetrates  through  holes  In  the  surface  and  Is  lost  Into  che  bulk  of 
the  material.  In  fact,  when  examined  In  chls  way,  the  Importance  of  surface  morphology  is 
most  striking  In  chat  the  same  relative  yields  are  also  found  for  the  low  Index  faces  of  Ag 
(18)  and  Au  (19).  The  placement  of  atoms  In  a  lattice  Is  clearly  an  Important  feature  thac 
Influences  che  yield.  This  notion  Is  also  consistent  with  Che  early  transparency  theories 
of  sputtering  (123)  which  could  be  qualitatively  used  to  explain  chs  relative  yields  for  che 
above  slcuadon. 


This  see  of  calculations  also  provided  considerable  insight  into  the  mechanism  of 
ejecclon  of  the  particles  from  each  face.  The  so  at  common  ejection  aechanlsa  found  in  these 
calculations  is  the  "up-dawn"  aechanlsa  where  the  prinary  ion  drives  the  target  acoa  down 
into  the  second  layer,  while  at  the  saaa  tiae  forcing  its  nearest-neighbor  to  eject.  The 
aooentua  can  be  transaitted  through  several  atoas  along  a  close-packed  row  by  this  process 
as  illustrated  in  Figure  14.  Thus,  it  is  clearly  sensitive  to  the  arrangeaent  of  surface 


Figure  34  -  Two  forms  of  surface  a  layer  sputtering  aechanlsa.  They  occur 
for  impact  points  which  lie  close  to  the  axis  of  a  dose- 
packed  chain  of  surface  atoms.  From  reference  93. 

For  all  three  target  orientations,  most  of  th<  ejected  atoms  are  found  to  originate 
from  the  surface  layer  as  jhovn  1..  Table  III  (142).  Note  that  from  (001),  virtually  every 
atom  originates  from  the  first  layer,  even  though  there  is  significant  movement  of  atoms  4 
layers  below  the  surface.  On  (110),  a  fairly  large  fraction  (~25Z)  of  the  atoms  are  ejected 
from  the  second  layer,  which,  because  of  the  open  structure  of  che  top  layer,  is  also 
exposed.  For  (111)  approximately  95Z  of  the  atoms  originate  from  the  first  layer.  These 
results  indicate  thac  Information  obtained  from  che  ejected  atoms  will  be  highly  surface 
specific. 

Table  III 

Origin  of  Ejected  Copper  Atoms  Resulting  From 
600  eV  Hormal  Incidence  Ar+  Ion  Impact 

Number  of  Ejected  Atoms 


< 


Origin  of  Ejected  Atom 
(Atomic  Layer  Number) 
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B.  Angular  Distributions  of  Atom  Elected  from  Clean  Single-Crystals 

A  second  structure-sensitive  factor  arises  froo  an  analysis  of  Che  angular 
distributions  of  the  secondary  particles.  As  Uehner  showed  aany  years  ago,  these 
dlscrlbucions  are  highly  anisotropic  and  very  clearly  reflect  the  surface  symmetry.  There 
have  been  aany  acceapts  to  explain  these  distribute  a.  One  such  explanation  la  chat  the 
ejection  occurs  along  dose-packed  lattice  directions  which  extend  deep  within  the  crystal 
(79).  This  Idea  nicely  explained  the  peaks  In  the  angular  distributions  but  required  chat 
chere  be  quite  a  bit  of  long  range  order  In  the  solid  even  during  Che  iapact  event.  As 
shown  In  Figure  1 ,  chat  requirement  seeas  a  bit  hard  to  swallow.  Although  controversy 
existed  concerning  these  "focusona"  for  aany  years,  the  molecular  dynamics  calculations 
of  Harrison  dearly  showed  chat  the  ejection  was  doalnaced  by  near  surface  collisons  rather 
than  those  from  beneath  the  surface  (76).  Similarly,  Lehmann  and  Slgnund  (143)  proposed 
that  the  spoc  patterns  could  be  foraed  by  the  synaetry  of  the  surface  layer  by  a  kind  of 
blocking  effect.  This  was  a  qualitative  model  and  was  never  extended  or  applied  In  a 
predictive  fashion. 

Recent  molecular  dynamics  studies  have  provided  Che  most  detailed  explanations  of  the 
Uehner  spots  and  have  suggested  a  number  of  new  experiments  alaed  at  elucidating  surface 
structure.  Using  che  large  micro-crystallite,  the  angular  distributions  have  been 
determined  with  sufficient  statistical  accuracy  Co  compare  directly  with  the  appropriate 
equipment.  These  distributions  are  shown  In  Figure  35a-c  Cu(OOl),  Cu(llO)  and  Cu  (111) 

(144)  where  each  ejected  atom  Is  plotted  on  a  flat-plate  collector  an  arblcrary  dlscance 
above  che  crystal.  For  chis  representation,  the  polar  angle  Is  determined  by  the  radial 
distance  from  the  cencer  of  che  plate  to  che  plotted  point  and  che  azimuthal  angle  Is 
determined  by  the  position  of  che  point  on  a  circle  drawn  about  the  origin  with  a  radius 
corresponding  Co  che  appropriate  polar  angle .  For  each  picture  che  complete  angular 
distributions  are  determined  after  the  Impact  zone  of  Irreducible  symmetry  (Figure  20)  Is 
unfolded  to  completely  encompass  che  target  atom  (77),  Note  chat  for  the  Cu(001)  surface, 
there  is  a  general  similarity  to  Uehner's  original  experiments  and  near  quantitative 
agreement  to  studies  by  Muscat  and  Smith  (145)  who  found  che  peak  In  the  polar  angle  spectra 
to  occur  at  44*.  The  nearly  six-fold  symmetry  found  on  the  (111)  face  and  the  diffuse 
rectangular  pattern  found  on  the  (110)  face  are  In  qualitative  agreement  with  experimental 
results  (1,18,23).  It  Is  still  dangerous,  however,  to  make  quantitative  comparisons  to  the 
aarly  experiments  for  the  reasons  already  discussed. 

Of  particular  (nearest  Is  Chat  che  feacures  In  the  calculated  angular  distributions 
change  rather  dramatically  when  only  chose  particles  with  rsther  high  kinetic  energy,  >20 
eV,  are  selected.  Aa  shown  In  Figure  3Sd-f  the  fraction  of  particles  ejected  along 
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(100)  (110)  (III) 


Figure  35  -  Angular  distributions  of  ejected  particles  for  the  (100),  (110), 

(111)  orlentatlone.  Each  ejected  atom  is  plotted  on  a  flat-plate 
collector  an  arbitrary  distance  above  the  crystal,  (a)-(c)  All 
ejected  particles;  (d)-(f)  only  chose  atoms  whose  kinetic  energy 
Is  greater  chan  20  eV.  The  numbers  on  che  ordinate  refer  to  Che 
polar  deflection  angle  given  in  degrees.  From  reference  144. 

preferred  crystallographic  directions  Is  considerably  enhanced.  The  lower  energy  particles 
tend  to  have  a  diffuse  angular  distribution  since  they  are  ejected  late  In  the  collision 
cascade  when  much  of  che  surface  order  is  no  longer  present.  This  is  a  general  observation 
and  one  which  can  be  applied  to  enhance  the  structure  sensitivity  of  any  angle-resolved 
measurement . 

Tracing  Individual  atomic  trajectories  yields  a  clear  picture  of  che  important 
scattering  mechanisms  In  che  high  kinetic  energy  regime.  First,  most  of  che  ejected 
particles  plocted  in  Figure  35d-f  arise  from  within  two  or  three  lattice  spacings  from  che 
Impact  point  and  suffer  only  a  few  scattering  events.  Second,  che  spacings  between  che 
surface  atoms  exert  a  strong  directional  effect  during  ejection.  The  fourfold  holes  on  the 
(001)  face,  for  example,  constrain  the  path  of  che  ejected  atoms,  and  their  trajectories 
proceed,  on  the  average.  In  che  (001)  plane  perpendicular  to  che  surface.  The  same 
mechanism  is  applicable  to  che  (111)  face;  che  three  elongated  lobes  arise  from  scattering 
through  the  threefold  holes  on  the  crystal  surface  that  do  not  have  second-layer  atoms 
directly  under  them.  In  ell  eeses,  the  atoms  ejected  from  the  alerocrystallite  originate 
from  the  top  layer.  Less  angle-resolved  structure  Is  observed  for  the  (110)  orlencaclon 
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since  che  rows  In  che  [100]  direction  are  sufficiently  far  aparc  co  allow  Che  atom  co  escape 
ac  virtually  any  angle. 

This  energy  selection  process  also  siaplifles  dramatically  the  number  of  operative 
ejection  mechanisms,  as  the  spot  Is  composed  nearly  entirely  from  Che  atoms  near  the  target 
atom  which  eject  frequently.  Since  these  particles  are  found  to  originate  near  the  target 
atom,  long-range  surface  order  should  not  be  required  co  produce  Che  spot.  A  further 
consequence  of  examining  only  che  higher  energy  particles  Is  that  che  trajectory  can  be 
stopped  after  a  shorter  number  of  time  steps  and  accurate  trajectories  can  be  obtained  by 
using  much  smaller  microcryscallltas. 

C.  Angular  Distributions  for  Adsorbate-Covered  Single  Crystals 

The  channeling  phenomena  observed  from  clean  surfaces  should  also  be  found  in  more 
complex  systems  such  as  for  metals  covered  with  a  chemisorbed  layer.  Only  a  few  attempts 
have  been  made,  however,  to  develop  a  comprehensive  theory  for  che  ejection  of  atoms  or 
molecules  from  reacted  surfaces.  The  specific  case  of  low-energy  He*  ion  bombardment  of 
nitrogen  adsorbed  on  tungsten  has  been  treated  using  the  3CA.  In  this  case  It  was 
postulated  that  che  nitrogen  atoms  could  only  be  desorbed  by  direct  collisons  with  the 
primary  Ion  (146).  Several  desorption  cross-sections  for  atoms  on  Ion  bombarded  surfaces 
have  also  been  calculated  using  che  BCA  via  che  computer  program  MORLAY  (147)  buc  agreement 
-  with  experiment  has  been  poor. 

Considerable  progress  In  quantitatively  describing  the  ejecclon  of  chemisorbed  atoms 
from  metals  has  been  made  using  molecular  dynamics  calculations.  The  first  model  system  Co 
be  explored  in  detail  was  oxygen  atomically  adsorbed  onto  Cu(001)  (99).  The  main  difficulty 
In  describing  any  situation  like  this  Is  co  develop  appropriate  Interaction  potentials  which 
describe  che  scattering  events.  Since  little  Is  known  about  chase  potentials,  early 
calculations  have  utilized  pair-wise  additive  potentials  for  adsorbates  which  have  the  same 
form  as  for  che  substrate,  buc  with  different  mass.  As  we  shall  see,  che  exacc  form  of  che 
potential  Is  noc  as  critical  as  the  atomic  placement  of  che  adsorbace  atom.  Thus,  In  che 
calculation,  che  geomecry  and  coverage  of  che  adsorbate  may  be  varied  over  a  wide  range  to 
test  how  these  quantities  Influence  ejection  mechanisms  and  ultimately  the  angular 
distributions. 

The  calculated  mechanisms  for  oxygen  ejection  In  this  case  are  quite  Interesting. 

There  are  two  common  pathways  for  oxygen  atom  ejection  which  Involve  direct  collisions  with 
the  primary  Ion.  The  Ar*  Ion  can  strike  the  target  oxygen  atom,  which  subsequently  reflects 
off  a  neighboring  Cu  atom  and  ejects.  Alternately,  che  Ar*  Ion  can  Initially  rcflecc  from 
the  target  copper  atom,  strike  a  nearby  oxygen  atom,  and  cause  it  to  eject.  Both  of  these 
processes  are  Illustrated  In  Figure  36a. 


MECHANISMS  FOR  ATOM  EJECTION 
-OXYGEN  ON  Cu(l 001- 
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Figure  36  -  Mechanisms  for  oxygen  atom  ejection  on  Cu(lOO).  From 
reference  99. 

Other  more  complicated  collision  sequences  can  be  Identified  as  Important  contributors 
to  oxygen  atom  ejection.  For  example.  In  Figure  36b,  the  copper  atom  adjacent  to  Che 
adsorbate  scarts  to  leave  the  solid,  but  finds  an  oxygen  atom  sticking  up  In  Its  pach. 
Depending  on  their  collision  angles  and  the  kinetic  energy  of  the  copper  atom,  both 
particles  may  or  may  not  eject.  An  oxygen  atom  can  also  be  desorbed  by  the  copper  acorn 
directly  beneath  It,  as  shown  In  Figure  36d,  although  it  Is  more  common  for  Che  copper  atom 
to  escape  the  solid  from  underneath  the  oxygen  without  ejecting  it,  as  shown  In  Figure  36c. 
The  reason  for  this  phenomenon  Is  that  only  a  small  percentage  of  the  energetic  copper  atoms 
eject  normally  to  the  surface. 

Since  che  lattice  surface  structure  dominates  Che  ejection  angle  of  Che  faster-moving 
particles  on  a  clean  surface,  It  appears  logical  chat  blocking  the  preferred  ejection 
directions  with  adsorbed  atoms  should  have  a  strong  effect  on  the  angular  distributions. 
Dslng  the  configurations  given  in  Figure  21,  for  example,  the  calculations  have  been 
performed  for  a  c(2x2)  overlayar  of  oxygen  on  Cu(001)  where  che  oxygen  Is  placed  In  an  acop. 
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2-rold  bridge  and  4-fold  bridge  bonding  geoaetry.  In  Figure  37  are  shown  the  relevant 
oxygen  distributions  for  the  atop  and  fourfold  bridge  site  (144,148).  As  It  turns  out,  the 
angular  distributions  for  the  Cu  substrate  stoas  arc  only  slightly  affected  by  Che  lighter 
oxygen  atoa  and  look  very  such  like  those  shown  in  Figure  35  for  clean  Cu.  For  oxygen  In 
Che  atop  site,  the  pattern  also  exhibits  ejection  angles  with  mainly  the  saae  symmetry  as 
Che  substrate.  For  Che  bridged  configuration,  however,  the  pattern  Is  rotated  45*  with 
respect  to  the  substrace  and  is  easily  distinguishable  from  Che  ocher  geoaetry. 
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Figure  37  -  Calculated  angular  distributions  of  oxygen  atoms  ejected  due  to 
Ar+  ion  bombardment  at  600  eV.  The  oxygen  originated  in  a 
c(2x2)  coverage  ou  Cu(001)  In  various  site  symmetries  and  heights 
above  Che  surface,  h.  The  nunbers  refer  to  the  polar  deflection 
angle,  (a)  A-top  or  linearly  bonded  site,  h  •  1.9  1. 

(b)  Fourfold  bridge  site,  h  ■  1 .2  i.  (c)  Fourfold  bridge  site, 
h  •  0.9  i.  From  reference  148. 

The  scattering  mechanisms  chat  give  rise  Co  these  angular  distributions  can  be 
discerned  on  an  atomic  level  from  an  analysis  of  the  trajectories.  It  Is  sufficient  to  note 
here  that  the  ejecting  substrate  atoms  most  strongly  Influence  the  directions  of  Che 
ejecting  adsorbate  atoms  and  chat  multiple  scattering  In  Che  overlayer  Is  of  secondary 
importance.  This  result  portends  chat  the  angular  distributions  of  the  oxygen  atom  might 
even  be  sensitive  to  Its  height  above  Che  surface  when  placed  in  a  given  geometry.  This 
hypothesis  la  borne  out  In  Figure  37c  where  the  oxygen  height  Is  reduced  from  1.2  1  to  0.9  1 
above  the  Cu  plane.  The  theory  suggest  thac  with  sufficiently  accurace  measurements  atomic 
positions  should  be  determinable  to  better  chan  tO.l  1. 

It  «>uld  be  delightful  If  experimental  methods  exlstad  to  teat  the  very  exciting 
predictions  regarding  these  angular  distributions.  With  Che  possible  exception  of  the  laser 
fluorescence  method,  however,  none  of  Che  methods  discussed  in  Section  II  have  anywhere  near 
the  sensitivity  to  detect  oxygen  atoms  from  an  oxygen  monolayer  before  destroying  the  sample 
Itself.  The  static  SIMS  method,  however,  provides  plenty  of  sensitivity  If  the  angular 
distribution  of  the  lone,  rather  chan  the  neutrals,  ere  measured.  This  approach  assumes 
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that  the  crejectori.ee  of  the  ions  and  neutrals  have  a  certain  correspondence  -  a  point  which 
appears  reasonably  valid  from  the  data  given  in  Section  III  B.  for  Che  C0/N1  system  015). 

Attempts  to  verify  the  predictions  of  the  molecular  dynamics  calculations  have  been 
carried  ouc  to  a  limited  extent  using  an  angle-resolved  SIMS  technique  ( 149) .  To  measure 
Che  angular  distributions  of  ejected  oxygen  and  copper  atoms,  a  cylindrical  shield  was 
placed  above  the  crystal  surface  as  shown  in  Figure  38.  The  polar  angle,  9,  for  ejected 
species  was  fixed  at  45*  with  two  appropriately  placed  apertures  while  the  azimuthal  angle, 
4,  could  be  varied  over  360*  by  rotation  of  the  crystal.  When  the  45*  electrostatic  sector 
in  front  of  Che  mass  spectrometer  was  utilized,  Che  kinetic  energies  of  the  ejected 
particles  could  be  roughly  selected  while  maintaining  a  large  enough  bandpass  necessary  to 
maintain  sensitivity.  This  experimental  configuration  will  produce  results  equivalent  to 
malting  a  circular  cut  of  che  spot  pattern  at  a  radius  corresponding  to  Che  45*  polar 
deflection  but  does  not  allow  the  entire  pattern  to  be  reproduced.  It  is  possible  to 
prepare  a  c(2x2)  overlayer  of  oxygen  on  Cu(OOl)  by  exposure  of  che  crystal  to  1200  L 
(1  L  ■  10“6  corr  s)  of  O2  at  25*C. 


CRYSTAL  I 

Figure  38  -  Schematic  representation  of  the  angle-resolved  SIMS  experiment. 

From  reference  149. 

The  observed  azimuthal  plots  for  Cu+  and  0"  ejection  from  these  samples  are  shorn  in 
Figure  39a.  The  fourfold  symmetry  of  the  (001)  orlentacion  is  clearly  evident  for  both 
species,  although  their  maxima  in  intensity  are  out  of  phase  by  45*.  This  result  is  only 
consistent  with  the  calculated  results  if  the  oxygen  has  adsorbed  in  the  fourfoLd  bridge 
site.  The  detailed  comparison  taken  using  a  fourfold  averaging  and  background  subtraction 
procedure  is  shown  in  Figure  39b,  assuming  che  oxygen  is  1.2  1  above  the  surface.  If  the 
calculation  had  been  performed  for  oxygen  placed  only  0.9  1  above  the  surface,  the  predicted 
curve  would  be  considerably  different  (114). 


59 


0  90  180 

^ ,  degrees 


4>,  degrees 


Figure  39  -  (a)  The  Cu+  and  0“  azimuthal  plots  as  recorded  from  Che  spectro¬ 

meter.  (b)  A  fourfold  average  of  the  data  In  (a)  with  the 
minimum  Intensity  subtracted  from  each  curve.  The  dashed  lines 
represent  the  calculated  curve  for  0  placed  In  a  fourfold 
hole  l .2  A  above  the  Cu  plane.  The  circles  represent  the  Cu+ 

Intensities  and  the  squares  represent  the  0"  intensities.  The 
primary  ion  is  1500  eV  Ar+  ac  a  total  dose  of  1-3  ions/cm^. 

From  reference  149. 

Of  crucial  importance  to  the  interpretation  of  these  spectra  Is  to  determine  how 
sensitive  the  calculated  patterns  are  to  the  vagarltles  of  the  input  potential  parameters. 
This  question  has  been  examined  in  some  detail  for  Che  O2  on  Ni(OOl)  system  by  Kapur  and 
Garrison  (97,150).  These  workers  found  chat  it  should  be  reasonably  safe  to  assign  the  type 
of  adsorption  site  (e.g.  atop  or  bridge-bonding)  even  if  the  potential  parameters  are 
uncertain  since  the  channeling  directions  are  basically  unaffected  by  changes  in  the  atomic 
sizes.  If  the  adsorbate  height  is  of  interest,  however,  they  found  chat  it  was  Important  to 
know  the  effective  size  of  Che  adsorbed  atom  at  the  scattering  energy.  Otherwise,  there 
could  be  a  large  uncertainty  in  the  determined  crystallographic  coordinates. 

Although  these  types  of  measurements  and  calculations  are  only  in  their  very  early 
stages,  the  preliminary  results  have  important  implications.  First,  the  experiment  Itself 
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provides  a  unique  situation  which  allows  detailed  comparison  between  experiment  and  theory. 
The  level  of  agreement  between  Ion  yield  measurements  and  calculated  neucral  yields  suggests 
that  their  Is  a  close  correspondence  between  the  trajectories  of  the  two  types  of  species. 
This  observation  has  Important  Implications  in  distinguishing  the  various  ionization 
mechanisms  proposed  In  Section  III.  And  finally.  If  details  of  the  experiment  and  the 
theory  can  be  further  refined,  the  approach  may  prove  to  be  a  powerful  new  method  for 
determining  the  positions  of  atoms  on  surfaces. 

V.  Cluscer  Formation  Processes 

A  second  critical  aspect  of  the  ion  bombardment  process  involves  the  fact  that  molecular 
cluscer  species  are  often  observed  to  be  ejected  from  the  surface.  As  discussed  in  section 
II,  a  variety  of  different  types  of  cluster  have  been  observed  from  a  wide  variety  of  sample 
subscraces.  For  adsorbates  on  clean  metals,  the  simplest  types  are  composed  of  Che  metal 
atoms  themselves  as  where  n  can  be  as  large  as  12  or  more  (46).  In  addition,  chere  can 
be  combinations  of  adsorbate  atoms  with  Che  metal  atoms  to  form  a  myriad  of  mixed  clusters. 
The  family  -  for  oxygen  covered  Si  where  n“0-4  and  s*0-3  represents  an  example  of  this 
type  (151).  In  addition,  there  have  been  many  cases  where  a  molecular  adsorbate  has  been 
observed  co  eject  intact,  as  is  the  case  for  fairly  large  organic  molecules  such  as  for 
p-aminobenzolc  acid  adsorbed  onto  silver  (152). 

There  are  many  fascinating  questions  regarding  the  formation  of  chese  clusters.  For 
example,  does  the  presence  of  a  given  molecular  fragment  indicate  that  species  actually 
existed  on  Che  surface?  Or  is  the  bombardment  process  too  catastrophic  to  preserve  the 
fragile  information  carried  through  chemical  bonds?  If  the  clusters  do  arise  from 
contiguous  surface  atoms,  can  their  presence  provide  key  information  regarding  che  local 
atomic  structure  of  surfaces  of  alloys  and  supported  metal  catalysts?  The  association  of 
NljCO’"  and  SiCO’'  co  bridged  and  linear  bound  CO  (Figure  15)  is  an  application  of  this 
concept.  Further,  is  there  any  point  co  measuring  che  angular  distributions  of  ejected 
clusters  and  can  they  be  used  to  probe  surface  structure?  And  finally,  how  does  che  charge 
on  the  cluster  lr  luence  its  formation  probability?  'dill  there  be  a  difference  in  cluster 
formation  mechanisms  for  positive  ions,  negative  ions  or  neutrals?  In  this  section,  we 
shall  probe  the  mechanism  of  cluscer  formation  and  examine  some  of  che  past  and  fucure  uses 
of  che  clusters  in  the  ion  bombardment  experiment. 

A.  Theoretical  Aspects  of  Molecular  Ejection 

Ever  since  the  initial  discovery  (44)  that  molecules  could  be  ejected  from  ion 
bombarded  solids,  the  genesis  of  cluscer  formation  has  been  the  subject  of  numerous 
investigations.  The  first  major  effort  arose  after  the  development  of  a  time-of-f light 
technique  aimed  at  measuring  the  energy  distribution  of  ejected  neutral  dimers  053) 
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provided  Che  energy  distribution  of  Kj  sputtered  fron  X.  With  this  data,  Konnen,  Tip  and 
DeVries  were  able  to  apply  the  results  of  the  transport  theory  to  obtain  a  theoretical  dimer 
energy  distribution  (48) .  By  assuming  that  the  atoms  In  the  dimer  arose  from  neighboring 
atoms  on  the  surface  and  by  utilizing  Thompson's  expression  for  the  energy  distribution  of 
the  monomer,  they  computed  the  normalized  probability  density  for  two  particles  which  have  a 
given  relative  energy  and  center  of  mass  velocity.  The  resulting  energy  distributions  were 
then  compared  co  the  experimental  data  for  K2  and  X  sputtered  from  a  polycrystalllne  K 
target  with  satisfactory  agreement  between  experiment  and  theory.  The  theory  has  also  been 
extended  to  include  the  formation  of  larger  clusters  (67)  such  as  those  observed  by 
Staudenmaler  (46).  The  success  of  their  predictions  led  these  workers  to  propose  that  the 
clusters  are  formed  by  recombination  of  Independently  sputtered  atoms  and  that  their  energy 
distributions  are  determined  by  the  £~Q  power  laws  of  the  monomer  distributions.  The 
concepts  have  never  been  fully  exploited,  probably  since  the  energy  distribution  formulae 
that  serve  as  the  basis  of  the  development  are  in  themselves  rather  crude  approximations  and 
are  not  easily  generalized. 

The  molecular  dynamics  calculations  have  been  extremely  successful  in  enhancing  our 
microscopic  understanding  of  how  the  clusters  form  and  in  helping  to  evolve  a  number  of 
geheral  considerations  In  developing  cluster  formation  theories.  The  first  study  regarding 
this  subject  was  completed  by  Harrison  and  Delaplain  (93)  who  investigated  the  bombardment 
of  Cu  single  crystals  by  600  aV  Ar+  using  the  classical  equations  of  motion  given  In  section 
II.  Although  chey  utilized  a  small  microcrystallite  and  were  only  able  to  sample  a  few 
Impact  points  on  the  surface,  using  Equation  22,  they  were  able  to  find  a  number  of  examples 
where  two  atoms  weie  above  the  surface  with  negative  total  energy.  This  condition  means,  of 
course,  that  In  the  absence  of  other  energy  transfer  processes  that  the  atoms  are  bound  onto 
a  stable  dimer.  Of  primary  Interest  In  this  study  was  that  the  atoms  that  formed  the  dimer 
were  not  atoms  that  were  contiguous  on  the  Cu  surface.  In  fact,  they  found  that  most  dimers 
form  from  next-nearest  neighbors  since  the  up-down  mechanism  (Figure  33)  reduced  the 
probability  that  neighboring  atoms  would  have  low  relative  kinetic  energies.  Thus,  the  study 
supported  the  work  of  Konnen,  Tip  and  DeVries  (48)  who  believed  that  the  dimers  formed  over 
the  surface,  but  questioned  their  assumption  about  the  proximity  of  atoms  that  form  a  dimer. 

Harrison's  early  landmark  calculations  were  mechanistically  very  important  but  they 
failed  to  provide  any  real  guides  for  the  experimentalists.  Mora  recently,  however,  the 
extent  of  the  calculations  have  been  sufficiently  expanded  and  refined  such  that  a  number  of 
additional  aspects  of  the  problem  can  be  examined  in  detail.  One  refinement  involves  the 
epproprlate  choice  of  interaction  potential  used  to  teat  for  dimer  formation.  The  solid 
pair-potentials  account  for  many  body  Interactions  but  do  not  have  the  correct  dissociation 


Holes.  For  example,  Cuj  hss  s  gas  phase  wall  depth  of  2.05  eV  (154)  whereas  the  bulk  pair 
potential  has  a  well  depth  of  only  0.48  eV  (93).  Since  the  tested  interactions  occur  above 
the  solid,  then  the  dlaer  pair  potential  should  be  oore  appropriate  to  use  In  determining 
cluster  stability  than  the  solid  pair  potential.  Other  Improvements  In  the  calculation 
relace  co  Che  Increase  In  computer  speed  available  from  Che  latest  computer  systems.  With 
this  technical  advance,  it  Is  feasible  to  generate  a  set  of  trajectories  with  statistically 
reliable  yields  for  many  types  of  clusters  using  mlcrocrystallites  that  are  large  enough  to 
contain  the  Ion  Impact  event  (155). 

tn  general,  the  calculations  Indicate  that  there  are  a  surprisingly  large  number  of 
clusters  that  would  be  expected  to  form  by  this  mechanism  as  has  been  indicated 
qualitatively  by  a  number  of  experiments.  Although  there  are  no  experiments  which  provide 
dimer  yields  from  single  crystals,  Oechsner  has  measured  the  Cu2/Cu  and  Guj/Cu  ratio  from 
polycrystalline  Cu  using  his  SNMS  technique  (29),  which  are  compared  to  calculated  values 
for  Cu(001)  in  Table  IV.  The  comparison  of  the  actual  yields  are  In  excellent  agreement  for 
Che  dimers  and  In  poor  agreement  for  the  crlmars.  The  calculated  number  could  be  too  high 
either  from  the  use  of  an  overly  strong  potential  to  describe  CU3  or  from  the  fact  that  some 
of  che  Cu3  molecules  could  dissociate  before  reaching  the  detector.  The  experimental  number 
might  be  somewhat  low  since  fraguentaton  of  CU3  during  Ionization  was  assumed  to  be 
negligible . 


Table  IV 

Theoretical  and  Experimental  Comparison  of  Cu2/Cu  and  CU3/CU  Ratio 
(001)*  Polycrystalline15 

Cui/Cu  .077  0.104 

Cujj/Cu  .014  0.001 

« 

•Calculations  performed  at  1000  eV  Ar+  Ion  Incident  at  9*0 *  using  a 
Moliere  potential  to  describe  the  Ar+-Cu  Interaction. 

'’Experimental  data  taken  for  1000  ev  Ar+  Ion  incident  at  9-0*  (29). 

There  has  also  been  found  to  be  a  strong  dependence  of  che  cluster  yields  on  surface 
morphology.  Recent  experimental  studies  have  been  completed  on  Nl(00l),  (110)  and  (111) 
where  che  relaclve  yields  of  Nl2+/Ni't‘  and  Ni3+/ Sf"  have  been  determined  to  be  In  the  same 
relative  order  as  the  monomer  yields  (156).  The  experiments  were  fortunately  completed 
under  conditions  appropriate  for  comparison  to  the  dynamics  calculation,  l.e.  ultrahlgh 
vacuum  with  low  Ar+  ion  doses  co  avoid  surface  damage.  Even  though  che  calculations 

ucllizad  attractive  potentials  for  CU2  rather  than  for  Nip  and  they  considered  only  neutral 

dimers  and  not  che  Ionic  species,  the  same  trends  In  the  calculations  were  observed  as  for 
che  experiment  (155).  Apparently  theae  other  effects  are  not  as  Important  In  Influencing 
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relative  cluster  yields  as  are  relative  positions  of  surface  atoms. 

These  Ideas  can  be  qualitatively  extended  to  explain  the  presence  of  even  larger 
clusters.  The  atomic  motion  that  leads  to  the  formation  of  a  pentamer  from  Cu(Hl)  has  been 
examined  In  detail  (94).  In  this  case,  as  shown  In  Figure  40,  there  are  three  collision 
sequences  Interacting  In  a  concerted  fashion  which  produce  5  atoms  moving  In  parallel  above 
the  surface.  Mote  that  the  atoms  that  comprise  the  Cuj  molecule  (atoms  #2, 3, 4,5  and  6)  do 
not  consist  of  contiguous  atoms  but  do  arise  from  a  highly  localized  region  of  the  crystal 


Figure  40  -  Schematic  representation  of  the  pentamer  formation.  Only  the 

atoms  Involved  In  the  ejection  of  the  Indicated  atoms  are  drawn. 

Other  atoms  are  located  at  the  Intersection  of  the  grid  lines  but 
are  not  shown.  The  actual  size  of  the  crystallite  used  in  the 
calculation  ts  considerably  larger  but  Is  not  shown  for  graphical 
clarity.  The  sizes  of  Che  acorns  are  arbitrary  to  allow  maximum 
viewing  through  the  crystal.  The  arrows  Indicate  the 
approximate  direction  of  motion  of  each  atom  during  the 
trajectory.  The  numbers  near  the  atoms  are  simply  used  as 
labels  (see  text).  From  reference  94. 

For  clean  metals,  Chen,  It  appears  chat  the  clutters  form  by  a  recombination  of  stoms 
over  the  surface  of  Che  crystal  while  the  atoms  are  still  within  Interaction  range  of  che 
solid.  Since  the  mechanisms  by  which  certain  correlated  collision  sequences  lead  to  cluster 
formation  appear  to  be  cryscal  structure  dependent.  It  makes  sense  to  examine  che  angular 
distributions  of  ejected  dusters.  The  fact  that  there  were  strong  angulsr  effects  was 
first  noted  by  Staudenmaler  (157)  for  W2+  ejected  from  ’><(110).  In  that  case,  he  found  chat 
the  Wj+  Intensity  maximized  at  similar  angles  as  as  che  V*  Intensity.  Those  experiments, 


however,  were  preformed  with  ISO  keV  Ar+  lone  under  heavy  dose  conditions  asking 
quantitative  coaparlsons  to  calculations  difficult. 


6  4 


A  aore  recent  study  on  Nl(001)  using  static  SIMS  with  2  keV  Ar+  Ions  produced 
qualitatively  slollar  results  ss  shown  In  Figure  41  <114).  The  azimuthal  distributions  for 
Ml*  and  Nlj"*-  peak  In  the  <100>  directions  (9*0*)  although  the  Nl2+  distribution  Is  more 
anisotropic  than  the  Ni*  distribution.  As  It  turns  out,  the  molecular  dynamics  calculations 


Figure  41  -  Experimental  angular  distributions  of  Nl+  and  Nlj4-  ejected 

from  Sl(OOl)  at  a  polar  angle  of  45*  i5* .  The  center-of-mass 
kinetic  energy  of  the  particles  is  between  10  and  50  eV.  Both 
curves  are  fourfold  averages  of  the  original  data.  The  Incident 
Ar+  Ion  has  2  keV  of  kinetic  energy  and  is  at  normal  Incidence. 
The  solid  is  at  room  temperature.  The  peak  counts  are  -900  and 
-500  counts/sec  for  the  Ni+  and  Nlj*  distributions,  respectively. 
The  <100>  azimuthal  directions  correspond  to  4*0* .  From 
reference  114. 


provide  very  almll  -  angular  distributions  for  SI  and  SI2  If  only  the  higher  kinetic  energy 
particles  are  considered.  The  results  of  over  1000  trajectories  are  shown  In  Figure  42. 

For  the  higher  energy  dimers,  the  calculations  predict  chat  at  9*0*  nearly  all  are  formed 
from  che  specific  mechanism  Illustrated  In  Figure  43.  With  this  process,  Chen,  the  observed 
dimers  should  originate  from  next-nearest  neighbors  along  the  close-packed  row. 

In  a  similar  vein.  It  Is  also  conceptually  possible  to  envision  certain  mechanisms  that 
dominate  che  ejection  process  when  the  Ar+  ion  is  Incident  on  the  crystal  ac  non-normal 
angles.  Foley  and  Garrison  have  predicted  that  the  next-nearest  neighbor  dimers  predominate 
If  the  Ar>  Ion  Is  Incident  at  9*0*  but  that  nearest  neighbor  dimers  become  more  Important  ac 
9*45*  and  4*45*.  This  effect  Is  shown  quantlcsclvely  In  Figure  44.  Similar  calculations  an 
a  two  dimensional  slab  of  SaCl  suggest  thet  It  Is  possible  to  peel  fairly  long  chains  of 
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Figure  42  -  Calculated  angular  distributions  of  N1  and  Hlj  at  a  polar  angle 
of  IS*  IS* .  The  center-of-mass  kinetic  energy  is  between  10  and 
50  eV.  Also  shown  Is  the  contribution  to  the  Nlj  distributions 
f roe  diners  formed  on  atoss  No.  1  and  No.  3  of  Figure  43.  The 
peak  Intensity  Is  110  and  12  particles  for  the  Ni  and  Nlj 
distributions,  respectively.  Froa  reference  114. 


contiguous  atons  from  the  surface  at  certain  Incident  angles  of  the  primary  Ion  (102). 

Thus,  although  rearrangement  clearly  occurs  as  the  atoms  that  fora  the  clusters  leave  the 
solid.  It  may  be  possible  to  evaluate  the  degree  of  rearrangement  by  selection  of  the 
appropriate  angles.  If  this  concept  were  extrapolated  to  alloy  surfaces  such  as  CuNi,  the 
relative  placement  of  the  alloy  components  should  be  able  to  be  determined  from  the  angular  ' 
distributions  of  Che  NI2,  Cuj  and  NICu  species. 

The  next  level  of  chemical  complexity  involves  cluster  formation  from  Ion  bombarded 
mmtal  surfaces  chat  have  reacted  with  simple  atoms  or  molecules  to  form  an  adsorbed  layer. 

Ne  examined  the  mechanisms  of  ejection  of  the  atoms  for  the  oxygen/Cu  situation  In  Section 
IV,  but  It  was  not  clear  whether  or  how  clusters  could  form.  We  now  consider  what  factors 
are  Important  In  Influencing  cluster  formation  for  two  different  situations;  chat  Is  oxygen 
adsorption  which  often  produces  atomic  oxygen,  and  carbon  monoxide  adsorpcion  where  Che  CO 
actaches  to  Che  surface  In  a  molecular  state. 

For  Che  cate  of  atomic  oxygen  adsorption,  we  consider  first  Che  possibility  of  forming 
Che  O2  molecule  via  Che  recombination  mechanism  proposed  to  exist  for  cluster  formation  from 
clean  surfaces.  It  is  possible  to  test  for  O2  formation  In  the  calculation  by  computing  all 
0-0  Interactions  between  ejected  oxygen  atoms  using  the  Interaction  potential  for  gas  phase 
O2  (99).  The  results  clearly  show  chat  Che  coverage  and  adsorption  site  syasecry  are  major 
factors  In  determining  which  clusters  are  likely  to  form.  The  most  spectacular  example 
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■«  *3  "  Hachanlsm  of  formation  of  th«  Nl2  dinar  which  preferentially 
ejecta  In  the  (100)  directions,  contributing  the  majority  of 
lntanalty  to  the  peak  in  the  angular  dlatrlbutlon.  (a)  Nl(OOl) 
showing  the  surface  arrangeaenta  of  atoaa.  The  numbers  are 
labels  while  the  X  denocas  the  Ar*  Ion  lapact  point  for  the 
aechanlsa  shown  In  figure  43b.  Atoaa  1  and  3  eject  as  Indicated 
by  che  arrows  forming  a  dinar,  which  la  preferentially  moving  In 
a  (100)  direction,  (b)  Three  dlaenslonal  rapressntatlon  of  a 
312  diner  format  lot.  procesa.  The  thin  grid  Unas  are  drawn 
between  the  naaraet-nelghbor  SI  atoaa  In  a  given  layer.  For 
graphical  clarity,  only  tha  atoaa  directly  Involved  In  the 
aechanlsa  are  shown.  Froa  reference  114. 
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field  of  dlaars  versus  original  separation  of  conatltutent  atoms 
on  che  surface.  The  nearest  neighbor  separation  on  copper  la 
2.56  1.  The  positions  of  tha  lines  are  the  separation  distances 
of  the  atoaa  In  the  surface,  (a)  9»4J*,  a-0*;  (b)  d-43* 

P-1 1 .23" ;  (e)  3-43*  ,  s-22.50*;  (d)  3-4  5*  ,  »-33.75*;  (a)  9.45* 
«"45* .  Froa  reference  92. 
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Involves  che  efface  of  coverage  on  che  O2  formation  probability.  For  a  p(2x2)  coverage 
(0.25  monolayer)  che  calculations  show  chat  the  oxygen  atoms  are  simply  coo  far  spare  Co 
have  a  significant  opportunity  to  collide  within  interaction  range  of  the  solid.  For  a 
c(2x2)  coverage  (0.3  monolayer),  however,  the  probability  of  Oj  formation  is  enhanced  by 
about  a  factor  of  4. 

This  effect  has  also  been  observed  experimentally  for  O2  adsorption  on  Ni(001)  (139). 

As  shown  in  Figure  43  Che  02~/0~  racio  increases  to  a  value  of  0.04  where  LEED  experiments 
indicate  the  p(2x2)  structure  predominates.  This  value  rises  to  a  maximum  of  0.16  at  the 
c(2x2)  coverage.  Thus,  the  02~/0~  ratio  increases  by  a  factor  of  4  as  the  coverage  is 
doubled,  similar  to  what  is  predicted  by  the  calculations  for  oxygen  in  a  bridge  position 
(138).  The  fact  that  che  0  atom  density  must  be  reasonably  large  to  observe  O2  emission  has 
been  used  in  che  analysis  of  SIMS  data  of  oxygen  on  Pb  (139).  In  this  study  both  the  O2”  to 
0-  ratio  and  che  Pb02~  to  PbO~  ratio  were  reasonably  constant  versus  oxygen  exposure.  'Even 
at  low  coverages  of  oxygen  atoms,  and  PbOo*  clusters  were  being  detected.  This  implies 
a  nucleaclon  mechanism  for  che  adsorption  of  the  oxygen.  If  the  oxygen  were  adsorbing 
randomly,  one  would  expect  che  Oj"  to  0"  ratio  to  significantly  increase  with  increasing 
oxygen  exposure. 


LANGMUIRS 

Figure  45  -  Experimental  Ion  yields  vs  oxygen  exposure  for  0~  and  che 
Oj'/O”  ratio  from  Nl(001).  From  reference  158. 

Many  other  types  of  molecular  dustsrs  have  been  observed  to  be  ejected  from  a  variety 
of  oxidized  metal  single  crystal  surfaces.  As  mentioned  previously,  oxidized  Mi  is  a  good 
example  with  Sl+,  Ni:+,  Ml3+,  01,  02*.  H10*,  Nl20+,  N102-  and  Si203“  all  observable  (151). 
The  calculaclonal  model  provides  insight  into  how  chess  dusters  srise.  For  the  case  of 
oxygen  on  Cu(001),  the  atomic  trajectories  clearly  show  that  all  types  of  clusters  form  over 
the  crystal  surfaca  after  the  components  of  the  cluster  eject  in  s  more  or  less  Independent 
manner.  The  mechanism  of  cluster  formation  is  completely  analogous  to  that  proposed  for  the 
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clean  metals.  Conalderable  rearrangement  of  atone  on  Che  surface  chac  fora  the  clueter  la 
possible;  Che  species  are  hoc  lifted  out  of  Che  surface  lncact,  although  considerable  local 
aconlc  order  is  preserved.  Since  the  constituent  atoms  of  a  cluster  originate  from  a 
localized  region  on  the  surface,  the  Incorporation  of  oxygen  Into  a  nickel  crystal  to  form 
N10{ 001)  should  affect  the  Ml  cluster  yields.  In  these  experiments  the  Ni3+  t rimer  yield  Is 
more  dramatically  reduced  chan  Is  the  Nlj*  dimer  yield,  since  the  exr'  '4  NIO  crystal 
statistically  has  few  Ni  atoms  In  the  localized  region  (151).  /Ill  of  •»  Nl-0  dusters 

observed  In  SIMS  have  been  fpund  in  Che  calculations  to  form  by  this  “over  the  surface' 
mechanism. 

No  matter  how  favorable  the  surface  morphology  Is  for  cluster  formation,  the  dusters 
will  not  be  observed  If  they  are  not  thermodynamically  stable.  One  measure  of  the  stability 
of  dimers  Is  the  well  depth  De  of  their  mutual  Interaction  potential.  The  yield  of  CuO 
diners  as  a  function  of  a  hypothetical  gas  phase  D,  Is  shown  In  Figure  16  ( 160) .  The  yield 


Figure  46  -  CuO  yield  vs  gas  phase  well  depth  De.  The  oxygen  atoms  were 

placed  In  a  2-fold  bridge  configuration  In  a  c(2x2)  coverage  with 
a  binding  energy  ■  0.75  eV.  A  total  of  110  trajectories  were 
calculated.  From  reference  160. 

increases  approximately  linearly  with  Da  for  Da  <  4  eV  and  then  plateaus.  The  leveling 
off  la  due  to  a  finite  number  of  slow  moving  atoms  thst  can  form  clusters.  If  the  intial 
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energy  of  Che  Incident  Ion  were  lncreaeed  so  thee  Che  local  number  of  acoaa  available  Co 
fora  clusters  were  Increased,  chen  che  linear  region  would  be  extended  to  larger  values  of 
D#.  Noes  chat  even  ac  De-0.1  eV  chare  are  two  pairs  of  acoaa  with  excreaely  saall  relative 
kinetic  energy  so  chat  two  bound  diaers  ara  observed. 

The  binding  energy  Ej,  of  an  adsorbaca  to  the  surface  also  Influences  the  yield  of 
cluscers.  The  effect  of  Ej,  on  the  CuO  dlaar  yield  Is  shown  In  Figure  47 .  In  all  cases  che 


Figure  47  -  CuO  yield  vs  binding  energy  E),.  Curve  *  Is  for  the  c(2x2) 

fourfold  bridge  configuration,  curve  b  Is  for  Che  e(2x2)  A-top 
configuration,  and  curve  £  Is  for  the  p(2x2)  A-top  configuration. 

A  well  depth  Da  of  2.7  eV~is  used  to  determine  che  yields.  Froa 
reference  160. 

qualitative  trend  Is  che  saae,  che  nuaber  of  CuO  diaers  decreases  with  Increasing  binding 
energy  of  che  oxygen  to  che  surface.  The  calculations  show  chat  che  0  atoa  yield  also 
decreases  with  Increasing  binding  energy  (99),  thus  che  CuO  yield  reflects  che  cotal  number 
of  oxygens  available  for  bonding.  For  all  che  calculations  shown  in  Figure  47,  the  Cu 
yield  la  nearly  constant,  tf  che  CuO  aolecules  were  ejecting  Intact,  one  would  expect  che 
opposite  crend  -  che  aore  strongly  bound  the  oxygen  la  to  the  copper  che  aore  CuO  diaers 
chat  should  be  found  ejecting. 

The  response  of  a  surface  to  Ion  boabardaanc  covered  with  a  oolecularly  adsorbed 
species  Is  aechanlsclcally  distinct  from  che  atomic  adsorbate  case.  For  CO  on  Nl(OOl),  for 
exaaple,  the  strong  C-0  bond  of  11.1  eV  and  the  weak  S1C0  bond  of  1.3  eV  help  co  keep  che 
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molecule  cogechar  during  ejection.  In  the  experimental  studies,  Che  ruin  peaks  Is  the  SIMS 
spectra  for  a  Mi(OOl)  surface  exposed  to  a  saturation  coverage  of  CO  are  Nl+,  Ni2+,  Nij+, 
NiCO*  Ml2CO+,  and  M^CO4-  as  shown  in  Figure  48*  All  Ions  show  a  snoot h  increase  in 
intensity  with  CO  adsorption  and  teach  saturation  after  2-1  CO  exposure  (0.5  <*onolayer 
coverage).  The  yields  of  C*\  0+.  SIC4,  and  StO+,  are  all  leas  than  0.01  of  the  Nl+  Intensity 


Figure  48  -  SIMS  spectrua  of  Vt(UI)  after  ZL  exposure  to  CO.  The  primary  Ion 
was  1000  eV  \t+  with  a  total  dose  of  10 l2  Ions. 

The  classical  dynamics  treatment  for  CO  on  Nl(OOI)  yields  results  which  are  In 
qualitative  agreement  with  these  findings.  Approximately  SOZ  of  the  CO  molecules  chat  ejec 
are  found  to  eject  Intact,  without  rearrangement.  The  formation  of  M1C0  and  M12C0  duscers 
have  been  observed  Co  form  over  the  surface  via  reactions  of  Ml  atoms  and  CO  molecules.  Mo 
evidence  has  been  found  for  MIC  or  M10  clusters  In  Che  calculations.  The  Ion  bombardment 
approach,  then.  Is  a  very  sensitive  probe  for  distinguishing  becwsen  molecular  and 
dissociative  adsorption  procssses. 

A  number  of  workers  have  attempted  to  identify  structural  relationships  found  using 
other  techniques  such  at  LEED  and  vibrational  spectroscopy  to  cluster  yields  In  SIMS.  The 
correlation  of  Ml2C0+  to  bridge  bonded  CO  and  NICO*  to  linear  bonded  CO  Is  an  example  of 
this  approach.  As  It  happens,  the  calculations  clearly  show  that  the  mechanism  of  cluster 
formation  Is  not  consistent  with  this  picture  sines  Che  clusters  form  over  the  surface  via 
atomic  collision*.  Furthermore,  recant  combined  LEED/ SIMS  results  Indicate  that  Che  cluste 
Ion  yields  ars  not  directly  ralacsd  to  the  adsorbate/  substrate  geometry  (161).  The  c(CxC) 
structure  of  CO  on  Nl(OOl)  with  all  the  molecules  In  the  atop  sit*  gave  the  same  Ml2CO'r.' 
.VICO4-  ratio  as  cha  compressed  hexagonal  DEED  structura  which  must  have  both  acop  and 
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bridge  bonded  CO  molecules.  It  Is  of  interest,  however,  that  Flelsch  et  al.  (106)  did  find 
a  change  In  the  MlCO^/Hi'*’  ratio  as  the  hexagonal  structure  began  to  fora.  Dawson  and  Taa 
(162)  found  that  the  secondary  Ion  yields  could  not  be  slaply  correlated  to  linear  or  bridge 
bound  species.  They  obtained  the  SIMS  spectra  for  c(2x2)-CO  on  Hi (001)  where  the  CO  Is  all 
linearly  bound  and  for  CO  on  a  Nl(001)  surface  presaturated  with  H2,  where  the  CO  is 
thought  to  be  exclusively  in  a  bridge  bound  state.  They  found  no  correlation  between  the 
presence  of  bridge  bonded  CO  and  the  formation  of  MlvC0+.  In  fact,  these  workers  concluded 
that  changes  In  Ion  yields  were  more  closely  related  to  Che  heat  of  adsorption  which  varies 
with  surface  coverage  than  to  any  surface  structural  effects. 

Although  most  of  the  above  arguments  are  qualitative  In  nacure,  they  do  serve  to  muddle 
the  notion  that  the  cluster  yields  reflect  Che  original  bonding  geometry  of  the  CO.  On  the 
other  hand.  It  Is  clear  chac  angular  distributions  for  atomic  adaorbases  are  very  sensitive 
to  che  surface  structure  so  It  Is  not  unreasonable  to  anticipate  similar  effects  for  the 
HI/CO  system.  Extensive  calculations  using  che  molecular  dynamics  procedure  (115)  have  been 
completed  for  the  atop  and  twofold  bridge  bonding  configurations  but  statistical 
considerations  have  restricted  the  analysis  to  only  the  HI  atoms.  As  shown  In  Figure  19, 
when  che  CO  Is  in  che  acoo  geometry,  the- calculated  Mi  distributions  peak  along  azimuthal 
directions  which  are  similar  to  the  dean  surface.  For  the  twofold  bridge  case,  however, 
che  CO  overlayer  tends  co  randomly  scacter  che  ejecting  HI  atoms  producing  a  much  different 
-pattern. 


Figure  19  -  Predicted  angular  distributions  of  Hi+  Ion  intensity  for  the  CO 
adsorbed  in  atop  and  twofold  bridge  sites.  The  numbers  on  che 
ordinate  refer  to  the  polar  deflection  angle  given  In  degrees. 
The  display  is  oriented  the  same  as  Figure  IS  so  chat  the 
horizontal  and  vertical  directions  correspond  to  .  The 
points  beeween  che  circles  contribute  to  the  azirauchal 
distribution  at  >45r7*.  Only  those  particles  with  kinetic 
energies  of  3U  eV  were  counted.  From  reference  115. 


There  has  bcsn  a  recent  attempt  co  confirm  this  prediction  using  the  new  angle-resolved 
SIMS  apparatus  discussed  earlier  (120).  The  experimental  results  for  the  azimuthal 
distributions  for  Si*  were  found  to  be  in  excellent  agreement  with  calculated  values  for  Si 
with  atop  bonded  CO  if  Che  calculated  values  were  corrected  for  the  presence  of  the  image 
force  (115,138).  Furthermore,  chare  was  poor  agreement  with  calculations  when  assuming 
bridge  bonded  CO.  Thus,  there  is  no  question  that  the  angular  distributions  are  sensitive 
to  the  surface  structure  of  the  adsorbed  overlayer. 

We  are  lefc  with  somewhat  of  a  dilemma  with  regard  Co  the  interpretation  of  Che  static 
SIMS  experiments  for  Che  CO  on  Sl(001)  case.  The  original  interpretation  of  the  changes  in 
the  Nl2C0+/SiC0+  ratio  by  Barber  et  al.  (53,54)  is  clearly  Incorrect  as  indicated  by  the 
mechanistic  Implications  of  the  molecular  dynamics  calculations.  In  addition,  ocher  workers 
have  been  unable  to  extend  the  predictions  of  their  ideas.  However,  it  is  also  clear  that 
the  scattering  angles  are  sensitive  to  surface  structure  so  there  may  well  be  reasons  co 
correlate  certain  types  of  molecular  clusters  to  certain  types  of  surface  structures, 
although  in  an  indirect  fashion.  These  generalizations  and  final  resolution  of  this  problem 
muse  still  await  further  experiments  and  calculations  on  this  very  challenging  system. 

It  is  of  Interest  to  analyze  che  internal  scate  distribution  of  the  ejected  CO 
molecules  (100).  Since  the  collisions  which  lead  to  ejection  are  often  energetic,  it  is 
most  appropriate  co  picture  che  molecule  as  dumbell  shaped,  rather  chan  spherical  or 
ellptlcal  as  would  be  che  case  for  very  low  energy  collisions.  Thus,  the  molecule  is 
generally  struck  on  one  end,  causing  it  Co  spin  rapidly.  Using  the  standard  techniques  of 
small  molecule  classical  scattering  (163),  the  available  internal  energy  can  be  partitioned 
between  vibrational  and  rotational  energy.  As  shown  in  Figure  50,  a  significant  fraction  of 
che  ejecting  CO  molecules  are  calculated  to  be  vlbratlonally  excited.  A  force-fit  of  che 
results  co  a  Boltzmann  distribution  produces  a  vibrational  temperature  of  3400-15000  K. 

Note  that  che  comparison  to  Boltzmann  behavior  is  rather  poor,  indicating  che  system  is  far 
from  thermodynamic  equilibrium.  The  calculation  also  shows  that  che  rotational  states  are 
highly  excited,  with  an  almost  random  population  of  states  to  J  values  as  high  as  100.  The 
effective  rotational  and  vibrational  temperatures,  then,  are  vastly  different. 

Similar  experimental  evidence  for  this  type  of  lncernal  stace  distribution  has  recently 
been  reported  by  Lineberger  and  co-workers  (164).  They  examined  the  fluorescent  spectra  of 
sputtered  Oj*  and  FeO"  and  found  chat  the  vibrational  temperature  was  approximately  5000  X 
in  both  cases.  They  further  noted  chst  the  electronic  temperature  and  che  vibrational 
temperature  were  quite  different,  again  arguing  against  a  thermal  equilibrium  model.  Their 
spectrum  is  shown  in  Figure  51. 


VIBRATIONAL  STATE  DISTRIBUTION 


Figure  50  -  Vibrational  state  distribution  for  the  molecular  ejected  CO 

species.  The  temperature  corresponding  to  the  line  of  steepest 
slope  Is  3450  K,  while  the  temperature  corresponding  to  the  other 
line  Is  15400  K.  From  reference  94. 
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Figure  51  -  03“  photoelectron  spectra  obtained  with  the  gas-discharge  Ion 

source  and  the  sputter  Ion  source.  The  sharp  break  In  the  Oj 
hoc  band  progression  at  v«4  occurs  because  this  state  la  no 
longer  stable  with  respect  to  Oj  *  e.  From  reference  164. 
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From  Che  varlecy  of  Che  above  studies,  1C  Is  possible  co  puc  forch  several  general 
mechanisms  of  cluster  formation  which  should  be  useful  In  the  Interpretation  of  mass 
spectra.  Flrsc,  for  systems  where  the  atoms  in  the  solid  have  no  preferred  molecular 
ldenclcv,  the  recombinacion  mechanism  muse  surely  predominate.  This  situation  would  exist 
for  pure  metals,  alloys  or  chemically  reacted  surfaces  where  the  adsorbate  Is  not  too 
strongly  bonded  co  che  substrate  and  interacts  simultaneously  with  many  substrate  atoms. 

From  che  variety  of  the  above  studies,  It  Is  possible  to  put  forth  several  general 
mechanisms  of  cluster  formation  which  should  be  useful  in  the  Interpretation  of  mass 
speccra.  First,  for  systems  where  the  acorns  in  che  solid  have  no  preferred  molecular 
Second,  If  che  solid  contains  molecular  entitles  with  strong  Intramolecular  forces  but  weak 
lncermolecular  forces.  It  Is  clearly  possible  that  clusters  may  eject  Incact  by  lacclce 
fragmentation.  Examples  of  this  type  would  Include  molecular  solids  like  Ice  (165)  and 
molecular  adsorbaces  on  metals  where  the  adsorption  bond  strengths  are  much  lower  than  Che 
bond  strengths  within  the  adsorbing  molecule.  Finally,  It  is  possible  that  clusters  can 
form  by  a  combination  of  both  of  the  above  mechanisms.  For  the  case  of  Ice,  for  example, 
each  Individual  water  molecule  may  eject  Intact,  but  then  form  water  multlmers  over  che 
surface  by  recombination. 

B.  Molecular  Clusters  and_  the  Ionisation  Problem 

as  It  has  been  applied  so  far,  the  molecular  dynamics  model  Ignores  any  effects  related 
to  electron  exchange  processs  and  ionization.  It  Is  reasonable  to  question  the  validity  of 
Che  predictions  of  this  model  'o  view  of  the  fact  chat  positive  or  negative  Ions  are 
detected  in  the  SIMS  experiment.  A  number  of  workers  have  addressed  this  question  at  least 
In  a  qualitative  fashion,  and  the  observations  help  considerably  to  understand  che  role  of 
ionization  In  cluster  formation  processes. 

One  of  the  first  comments  on  chis  question  was  made  In  regard  to  differences  in  che 
interaction  potential  expected  for  charged  species  vs.  neutral  species.  For  oxygen  adsorbed 
on  Cu(OOl),  for  example,  the  molecular  dynamics  model  indicates  chat  the  0;  molecule  would 
be  expected  co  form  by  the  recombination  mechanism.  For  this  specific  system.  It  is  also 
possible  co  examine  che  effect  of  charge  state  on  che  dimer  formation  probability  since 
Morse  potential  parameters  are  known  for  Oj'*’  and  O2”.  By  putting  chese  parameters  as  given 
in  Table  V  into  Equaclon  29  and  using  Che  criteria  for  dimer  formation  In  Equation  22.  che 
probability  of  forming  a  dimer  between  an  0"  or  C&  ion  and  an  0  atom  can  be  deduced  (99).  Of 
particular  interest  is  that  the  number  of  calculated  dimers  differs  by  less  chan  a  factor  of 
2  for  che  three  cases  02+,  O2  and  O2".  The  most  striking  feature  of  the  results,  however. 

Is  chac  che  thermodynamically  most  stable  species,  02+,  produces  the  lowest  number  of 
dimers.  The  reason  for  this  feet  Is  thsc  the  range  of  the  potential  Is  the  smallest  for 
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Table  V 

Potential  Parameters  for  O2 ,  02 *  and  O2- 


Species 

De  («V> 

S  (i_1) 

Re 

°2 

5.12 

2.67 

1.21 

o2* 

6.48 

2.87 

1.12 

°2~ 

4.07 

2.11 

1.31 

thls  species.  Apparently,  Che  number  of  possible  two-body  Interactions  that  a  particular 
atom  might  experience  is  more  important  than  the  ultimate  stability  of  a  dimer  in  Its 
equilibrium  configuration.  If  It  Is  possible  to  generalize  this  Idea  further,  we  would 
expect  oppositely  charged  species  such  as  Ha'*'  and  Cl-  to  have  a  very  high  probability  of 
forming  an  NaCI  molecule  due  to  che  Infinite  range  associated  with  purely  electrostatic 
1 nee race Ions. 

A  test  of  the  formation  of  eluscer  ions  by  an  atom— Ion  collision  mechanism  was 
attempted  using  the  5(1(001)  surface  during  exposure  to  oxygen  (158).  In  this  case,  the  ions 

SI*,  NIO1,  0 1  and  Oj *  could  be  monitored  as  a  function  of  the  oxygen  dose.  The  results  as 

shown  In  Figure  52  are  striking,  as  the  shapes  of  the  ion  yield  vs.  coverage  curves  are 
quite  different  depending  upon  che  Ion  Involved.  Mote,  however,  chat  the  NiO*  yield  has  a 
similar  shape  as  che  Mi"*"  yield  whereas  che  S10”  and  0“  and  the  0n~  and  0*  curves  also 
seem  to  track  each  other.  If  one  Invokes  a  eluscer  formation  mechanism  aa  discussed  r.w.e, 
Chen  che  cluster  Iona  could  form  aa  a  result  of  collisions  between  neutral  and  Ionized 
atoms.  For  example,  che  NIO*  ion  should  form  by  Interactions  between  SI*  Ions  and  0  atoms. 
The  reaction  of  Ml  atoms  and  O'*  ions  would  not  be  important  due  to  the  low  levels  of  0* 

Iocs.  Furthermore,  If  the  number  of  ions  above  che  surface  is  much  less  than  che  number  of 
neutrals  (which  Is  che  case  for  most  materials  txcepc  perhaps  alkali  halides),  che  Intensity 

of  the  resulting  eluscer  Ion  will  be  controlled  by  che  amount  of  the  corresponding  atomic 

Ion.  Thus  for  SI,  che  clusters  would  form  by  the  following  reaction: 


Si*  + 

0 

- > 

S10* 

(64) 

surface 

Si  + 

0 ■ 

- ► 

N10* 

(65) 

surface 

0*  + 

0 

- » 

On* 

(66) 

surface 

The  mechanism  of  eluscer  formation  for  ions  would  therefore  be  very  similar  to  chat  of  the 
neutrals . 

This  Idea  was  examined  in  more  detail  by  Yu  (109)  in  a  sec  of  elegant  experiments 
• '  .  .  c 1  - g  le-osclon  on  oxidized  T1  surfaces.  Before  Cs  deposition,  he  observed  using 
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Figure  52  -  Experimental  Ion  yields  vs  oxygen  exposure  for  Nl(QOl)  bombarded 
by  2-keV  Ar*  Ions.  The  ordinate  gives  Che  number  of  lons/s 
detected.  From  reference  158. 

static  SIMS,  the  Ions  Tl+,  710*,  7102*  and  TIOj*.  If,  for  example,  the  710*  species  formed 
by  recombination  as 

71+  +  0  - -»  Tiff*-  (67) 

surface 

then  suppression  of  the  71'*'  signal  by  Cs  deposition  should  suppress  Che  710*  signal  In  a 
similar  fashion.  This  experiment  Is  superior  to  Che  0;  on  N1  situation  since  only  small 
quantities  of  Cs  (<0.25  monolayer)  need  to  be  deposited  to  change  the  Ion  yield  by  nearly 
three  orders  of  magnitude.  In  addition,  the  oxygen  coverage  could  be  kept  fixed  during  the 
enclre  analysis.  His  results  showed  that  for  the  positive  Ions  710*,  TiOj*  and  7103*  that 
the  recombination  model  was  clearly  valid.  For  the  negative  Ions,  however,  he  found  strong 
deviations  from  the  predicted  behevior.  Because  of  this  deviation,  he  concluded  that  these 
species  probably  formed  directly  by  lattice  fragmentation.  It  la  of  Interest  that  he  did 
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noc  consider  reactions  of  Che  cype 


Ti+  +  0*  +  0-  *  T102"  (68) 

or  Tl+  +  0  +  0"  +  0“  *  TiOj"  (69) 

In  his  analysis.  These  pathways  may  indeed  be  favorable  due  to  the  long  range  attractive 
Interaction  char  exists  between  oppositely  charged  atoms. 

Yu  has  proposed,  however,  that  since  Che  electron  affinity  of  negatively  charged  ions 
is  usually  low  that  unless  Che  duster  constituents  are  in  ground  vibrational  states  that 
electron  detachment  is  very  favorable  (166).  To  produce  ground  state  molecules  he  feels  it 
is  necessary  to  Invoke  Che  condition  that  they  form  from  nearest  neighbors  on  the  surface. 
Thus  for  the  Ti-oxygen  system,  Yu  proposes  that  the  positive  ions  form  by  recombination 
while  the  negative  ions  are  produced  direccly  by  lattice  fragmentation. 

Another  possible  mechanism  for  cluster  ion  formation  has  been  suggested  by  Snowden  and 
co-workers  (167).  In  this  scheme,  which  they  tested  by  examining  S2+  ejection  from  an  ion 
bombarded  nitrogen  implanted  silicon  target,  the  electron  is  assumed  to  be  emitted  after 
formation  of  the  N2  molecule  aa 

S  +  S  - 9  S2  - -f  »2+  +  le*  (70) 

surface 

Aa  shown  in  Figure  S3,  che  first  step  (I)  requires  the  ejection  of  two  nitrogen  atoms  which 
form  an  S2  molecule  by  the  recombination  mechanism.  In  the  second  step  (II),  the  excited  Sj 
molecule  vibrating  within  its  molecular  potential  may  cross  to  the  No*  potential  curve  by 
emitting  an  electron.  This  process  can  be  detected  in  step  (III)  since  the  excited  N2+ 
molecule  decays  to  che  ground  scats  via  a  measurable  photon.  They  provided  evidence  for 
this  mechanism  by  determining  chat  specific  rocacional  levels  of  N2*  uere  populated  during 
step  (II).  The  proposed  mechanism  has  considerable  meric  since  essentially  unbound  So 
molecules  can  be  stabilized  as  by  electron  ejection.  It  has  so  far  only  been 
demonstrated  to  occur  in  this  one  Instance,  and  it  remains  to  be  seen  how  common  such 
processes  really  are.  Further,  the  fact  that  they  also  found  a  continuum  of  vibrational  and 
rocatlonal  states  that  uere  populated  Implies  that  a  fraction  of  the  So*  molecules  may  form 
by  ion-atom  collisions. 

Thus  it  appears  that  the  relationship  between  ionization  and  cluster  formation 
introduces  s  few  additional  twists  whan  compared  to  cluster  formation  of  che  neutrals.  Many 
of  the  basic  ideas  generated  by  the  trejectory  calculations,  however,  still  appear  to  be 
valid  and  continue  to  provide  valuable  Inelghc  inco  che  cluster  Ion  formation  process. 


—  , 
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Figure  S3  -  Schematic  Illustration  of  the  proposed  sputtered  molecule 

formation  mechanism.  The  first  step  (I)  requires  the  sputtering 
of  two  nitrogen  atoms  approaching  on  an  appropriate  collision 
trajectory.  The  second  step  (II)  requires  a  molecular  potential 
curve  crossing  followed  by  electron  loss  within  a  few  molecular 
vibrations  (the  lifetime  of  the  82  state).  In  che  third  step 
(III)  the  MV*-  molecule  decays  to  the  ground  state  via  photon 
emission.  If  the  electron  loss  in  II  does  not  occur,  che  82 
molecule  generally  predissociates  and  no  is  formed.  From 
reference  5,  p.  298. 


VI.  Molecular  SIMS 

The  detection  of  molecular  ions  by  static  SIMS  which  are  characteristic  of  the  original 
sample  has  important  implications  in  chemical  analysis.  The  approach  offers  a  complementary  i 

mass  spectrometric  method  for  vapoi 'zing  and  ionizing  nonvolatile  or  thermally  unstable 

1 

compounds.  For  example,  in  traditional  mass  spectrometry,  a  molecule  is  usually  volatilized  j 

by  heating,  and  ionized  by  electron  Impact.  For  molecules  like  amino  acids,  however, 

heating  usually  produces  8H3  and  CO2  whereas  ion  bombardment  may  directly  produce  the 

desired  molecular  ion.  In  ocher  words,  che  sub-picosecond  time  scale  of  che  ion  impact 

event  is  sufficiently  fast  to  beat  out  thermal  rearrangements  chat  occur  on  the  nanosecond 

time  scale. 

The  analysis  of  molecular  solids  by  SIMS  daces  to  the  early  investigations  of 
Bennlnghoven  on  metal  surfaces  covered  with  rather  ill-defined  organic  overlayers  (168).  * 

From  chase  observations,  he  and  his  co-workers  found  that  even  for  very  delicate  organic  j 

molecules  such  as  amino  adds,  ion  bombarded  powders  of  the  sample  could  produce  very  I 

well-defined  spectra  (51).  For  example,  in  Figure  54  is  shown  both  the  positive  and  1 

_  --  -  -  _  ---  -  -  ,  | 
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Figure  54  -  Original  positive  and  negative  secondary-ion  spectrum  of  cysteine 
on  Ag.  (Primary  Ion  current:  2’ 10"^  A.  bombardment  target  area: 

0.1  cm2.  From  reference  51. 

negative  Ion  spectra  of  cysteine  deposited  onto  a  sliver  foil.  The  main  features  of  the 
spectra  are  the  presence  of  the  molecular  Ion  peaks  at  (M4-l)+  and  at  (M-l)~  and  the  rather 
simple  fragmentation  pattern  associated  with  the  lower  mass  ranges.  The  ultimate 
sensitivity  of  the  method  Is  <10”12  grams,  making  SIMS  a  potentially  very  sensitive 
analytical  tool  (169,170).  The  cysteine  example  is  a  particularly  powerful  one  since  it  is 
chermally  unstable  and  the  mass  spectroscopy  of  amino  acids  is  very  tricky. 

Many  ocher,  experimental  studies  have  now  boon  reported  which  indicate  that  SIMS 
provides  interesting  information  from  almost  any  molecular  solid.  It  is  not  our  intent  to 
review  this  field  as  it  has  recently  been  done  quite  capably  by  others  (171,172).  There  are 
a  number  of  mechanistic  features  of  these  experiments  which  are,  however  related  to  many  of 
the  ide.a  presented  in  earlier  tactions  upon  which  we  will  focus. 
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A.  A  Modal  Study  with  Benzene 

A  good  starting  point  in  crying  to  understand  Che  ejection  of  large  aoleculea  la  to 
begin  with  an  appropriate  aodel  system  -  In  our  case  benzene.  This  aolecule  Is  simple  enough 
to  be  examined  theoretically  yet  large  enough  to  be  comparable  In  size  to  many  of  Che  more 
complex  systems  currently  of  interest  to  chemists.  The  static  SIMS  spectrum  of  this 
compound  has  been  examined  under  a  variety  of  conditions.  Mlchl  and  co-workers  first 
reported  the  spectrum  of  solid  benzene  at  low  mass  (173)  and  was  quickly  extended  to  higher 
mass  by  Rabalals  ec  al.  (165).  The  spectrum  la  characterized  by  a  series  of  hydrocarbon 
dusters  of  the  stoichiometry  CnHx  where  n«l  to  30  and  x  is  on  the  order  of  n  for  nClO  and 
on  the  order  of  n/2  for  n>10.  The  spectrum  Is  dearly  dominated  by  recombination  of 
ejecting  cluster  sub-units,  making  direct  identification  of  benzene  itself  nearly 
impossible.  For  example,  Che  following  types  of  gas  phase  reactions  have  been  proposed  co 
explain  Che  formation  of  the  larger  carbon  fragments  (165): 


C3H34,  +■  *  C9h7+  +  H2 

(71) 

C5H3+  +  C6H6  -  CnH7+  ♦  H2 

(72) 

c5h5+  +  *  cllH9+  +  a2 

(73) 

A  dramatically  different  spectrum  is  obtained  by  allowing  benzene  co  chemisorb  on  a  Ml(OOl) 
surface  (174,175).  For  this  situation,  it  is  believed  chat  Che  molecule  forms  an  ordered 
c(4x4)  (1/8  monolayer)  overlayer  with  the  benzene  aolecule  laying  flat  on  the  surface  in  a 
dose-packed  arrangement  (176,177).  The  geometry  is  sketched  schematically  in  Figure  55. 

The  SIMS  spectrum  consists  almost  entirely  of  the  molecular  ion  complexed  with  a  Ml  atom  as 
NIC^H^'4’  as  shown  in  Figure  56.  Only  weak  fragmentation  patterns  are  observed.  Apparently, 
the  benzene  molecules  are  far  enough  apart  such  chat  recombination  of  the  ejected  molecule 
occurs  maidy  with  Ni+. 

One  may  quesclon  how  these  observations  can  be  rationalized  with  Che  mechanistic  ideas 
discussed  In  Section  V.  How,  in  the  case  of  benzene,  can  12  or  more  atoms  be  ejected  from  a 
solid  retaining  molecular  structure  information  whan  bombardment  energies  are  many  times 
greater  than  chemical  bond  strengths?  This  situation  has  recently  been  modeled  using  the 
molecular  dynamics  treatment  to  see  if  Che  theory  could  be  used  to  rationalize  Che 
experimental  result  using  the  structure  shown  in  Figure  55  (101)-.  The  model  was  set-up  by 
m^z  a  binding  energy  of  1.7  eV  between  Che  entire  benzene  molecule  and  the  Mi  surface. 

The  majority  of  Che 
rather  Chan  into  the 
ignored.  With  this 
with  smaller  amounts 


tlon  strjjpfth  was  arbitrarily  put  into  Che  Ml-C  interaction 
-H  incerriVlon.  Also,  Interaction  between  benzene  molecules  was 
approach,  W  analysis  indicated  very  high  probability  of  C$H$  ejection 
of  fragmm 
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also  found  where  the  species 


Its  of  the  type  CgH;,  CjHj.CjH?  and  CH.  One  case  of  MIC^K^  was 
is  formed  by  recombination  of  a  Ml  atom  with  a  unit. 


Nl(OOl)  with  a  c(4x4)  ovarlayer  of  benzene,  .  The  dashed 
triangle  is  the  inpact  zone  for  nornal  incidence  ion  bombardment. 
The  numbered  aeons  correspond  to  the  atons  shown  in  Figure  57. 

The  x  is  the  inpact  point  for  the  Ar*  ion  which  leads  to  the 
motion  depicted  in  Figure  57.  The  circle  around  each  C6 Rg  is  the 
radial  extent  of  che  hydrogen  position,  2.5  K.  From  reference  101 
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SIMS  spectrum  of  Hid  11)  exposed  to  31.  of  benzene  at  300  K 
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The  calculations  indicate  that  there  are  three  factors  that  favor  ejection  of  molecular 
fragments  (101).  First,  a  large  molecule  has  many  Internal  degrees  of  freedom  and  can 
absorb  energy  from  an  energetic  collision  without  dissociating.  Second,  the  more  massive 
framework  of  a  large  organic  molecule  can  be  small  in  size  compared  to  a  metal  atom,  thus  it 
is  possible  to  strike  several  parts  of  the  molecule  in  a  concerted  manner  so  that  the  entire 
molecule  moves  in  one  direction.  Finally,  by  the  time  the  organic  molecule  is  struck,  the 
energy  of  che  primary  ion  has  been  dissipated  so  chat  the  kinetic  energies  are  tens  of  eVs 
rather  chan  hundreds  or  thousands  of  eVs.  k  typical  ejection  sequence  for  benzene  Is  shown 
in  Figure  57 . 


0  fs 


Figure  57  •  Ejection  of  a  molecule  by  che  simultaneous  collision  with 

two  Ml  atoms.  Only  che  species  (one  Ar+  ion,  3  Hi  atoms  and  one 
molecule)  directly  involved  in  che  collision  are  shown.  The 
grid  lines  are  drawn  between  nearest  neighbors  in  each  plane, 
thus  a  Si  atom  is  initially  situated  ac  each  intersection  of 
lines.  The  sizes  of  the  atoms  are  arbitrary.  The  elapsed  else 
during  che  collision  is  shown  in  the  lower  right  corner  in 
femtoseconds  (10*13  s). 

e)  Initial  positions  of  the  atoms.  The  benzene  is  being 
viewed  from  che  side. 

b)  Positions  as  the  2  Si  atoms  are  about  to  collide  with 
benzene . 

c)  Distortion  of  the  ejected  benzene.  From  reference 
101 . 
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It  should  be  noted  here  that  the  SIMS  spectrum  of  water  exhibits  similar  variations  In 
cluster  type  with  changes  In  Che  sample  matrix.  The  spectrum  of  Ice,  for  example,  consists 
calmly  of  the  series  of  ions  CH^O^rC1"  with  n  up  to  at  least  31  (165).  When  wacer  is 
adsorbed  oa  a  cooled  meeal  substrate,  the  M(H20)o*  Ions  dominace  che  spectrum  (178). 

8.  Matrix  Effects 

The  comparisons  of  che  SIMS  spectra  of  benzene  under  the  different  matrix  conditions 
(solid  benzene  and  benzene  adsorbed  on  Ni(OOl))  suggest  that  fragmentation  and  recombination 
processes  can  be  strongly  influenced  by  sample  preparation.  This  fact  has  been  further 
demonstrated  by  a  number  of  workers.  Jonkman  and  Mlchl,  for  example,  utilized  a  liquid 
He  crystal  to  trap  che  organic  molecule  In  solid  argon  in  an  attempt  to  simplify  che  spectra 
(179-181).  For  solid  CH4,  they  found  hydrocarbon  fragments  to  masses  assigned  to  CqHyr  and 
higher.  Uheu  CH4  was  diluted  Into  an  argon  matrix  {1:5001  at  15  K,  CH4+  was  che  highest 
molecular  weight  species  observed  (180).  These  effects  are  Illustrated  In  Figure  58. 


(b) 


Flyura  38  -  (a)  SIMS  spectrum  of  neat  methane  at  15  K  using  l  keV  He*  tons. 
’’  (b)  SIMS  spectrum  of  methane  diluted  In  an  argon  matrix  {1:5001 

at  15  K  using  1  keV  He*  ions.  From  reference  130. 
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From  reaction*  of  the  type  given  in  Equation*  71-73,  it  1*  obvious  that  the  neutral 
aolecule  can  be  converted  into  an  ion  by  recombination  reaction*.  Cook*  and  co-workers  have 
taken  advantage  of  this  idea  by  mixing  large  quantities  of  metal  ions  such  a*  Ll+  and  K* 
into  the  sample  which  results  in  intense  eatlonlzed  molecular  ion  peaks  (182,183).  These 
workers  also  found  chat  p-aminobenzolc  acid  prepared  as  a  film  on  Ag  foil  formed  the 
argencaced  molecular  ion  (Ag+M)+  (182)  as  was  found  for  benzene  on  Ml(001). 

C.  Alkali  Halides 

The  situation  becomes  more  complex,  but  nonetheless  fascinating,  when  the  particles 
forming  the  cluster  possess  long  range  attractive  forces.  For  Che  alkali  halides  MA, 
clusters  of  the  type  (MA)nM+  and  (MA)nA~  are  the  most  commonly  observed  species,  with  n 
increasing  co  numbers  greater  than  20  (178,  184-186).  This  situation  is  illustrated  in 
Figure  39  for  the  positive  ion  spectrum  of  Csl  (184).  Presumably,  these  clusters  form  by 
recombination  of  an  If*  and  A"  over  the  surface,  a  very  favorable  process  due  co  the  infinite 
range  of  the  interaction.  The  (HA)  molecule  can  then  be  sequentially  eatlonlzed  or 
anlonized  over  che  surface  of  the  crystal  co  build  up  units  contained  in  the  observed 
cluster  ion.  The  energy  distributions  tend  to  be  narrower  than  the  monomer  ions  (187). 
Further,  there  are  certain  cluster  sizes  which  exhibit  unusual  stability.  The  peak  at 
n-13 ,  for  example,  is  thought  co  be  associated  with  a  rhombohedral  structure  containing 
27  atoms  arranged  in  a  3x3x3  structure. 

0.  Surface  Reactions 

In  general,  che  SIMS  spectra  of  monolayers  adsorbed  on  metals  appear  to  produce 
dusters  which  often  minimize  complexities  introduced  by  racomblnacloo  processes.  This  fact 
has  inspired  considerable  research  aimed  toward  utilizing  SIMS  co  actually  monitor  surface 
reactions  that  may  be  occurring  on  surfaces.  Bennlnghoven  and  co-workers  have  emphasized 
the  sensitivity  of  SIMS  co  hydrogen  by  studying  the  reaction  of  Hj  with  02,  C2H4  and  C1H2  - 
polycryscalline  HI  (188-190).  For  tha  Hj  +  reaction, for  example,  they  were  able  to 
detect  an  0H~  species  at  room  temperature,  a  result  similar  co  chat  found  for  che  adsorption 
of  H2O  on  Ml  under  similar  conditions  (189).  Other  workers  Have  followed  che 
dehydrogenation  of  ethylene  at  120  K  on  Mi(lll)  (192).  In  a  r«re  complicated  example,  it 
has  recently  been  observed  chat  thiophene  adsorbed  on  Ag  undergoes  a  self  hydrogenation  to 
form  cetrahydrothiophene  (193). 

*•  Molecular  SIMS  and  Related  Methods 

The  search  for  new  methods  of  introducing  ions  from  non-volatile  samples  into  a  mass 
spectrometer  has  been  of  interest  for  many  veers.  There  is  no  question  that  static  SIMS  has 
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Figure  59  -  Secondary  Ion  mass  spectra  of  Csl . 

Figure  courtesy  of  J.E.  Camp ana. 


found  a  niche  In  this  field  In  the  analysis  of  biochemical  materials.  Inorganic  and  organic 
solids  as  well  as  salts  and  metals.  It  Is  Interesting  that  there  are  currently  two  other 
approaches  of  very  similar  fundamental  origin  under  current  development.  Macfarlane  and 
co-workers  have  discovered  that  by  bombardmant  of  a  thin  mylar  foil  from  the  back  by  the 
fission  products  of  a  252cj  nucleus  that  molecular  Ions  are  ejected  from  the  front  of  the 
foil  (194,195).  This  technique  has  been  referred  to  as  pxasma  desorption  mass  spectrometry 
and  Is  compared  schematically  to  SIMS  In  Figure  60.  We  believe  that  the  same  general 
mechanisms  of  cluster  formation  are  applicable  to  this  method  as  well.  Although  the 
bombardment  energy  is  in  the  MeV  range  and  energy  loss  processes  are  almost  certainly 
electronic  In  nature,  any  subsequent  collision  cascade  emlnatlng  out  from  Che  fission  crack 
would  lead  to  the  ejection  of  clusters  In  a  fashion  nearly  identical  eo  SIMS.  The  fact  chat 
the  (MH)'**  and  (M-l)“  ions  dominate  Che  PDKS  speccrum  Is  qualitative  evidence  In  support  of 
this  view  (194).  of  further  Inceresc  Is  chac  a  cime-of-flighc  mass  spectrometer  has  been 
developed  co  detect  Che  Ions.  The  spectrometer  picks  up  Uj  timing  pulse  from  Che 
Individual  252^/  fission  event.  This  feature  allows  very  high  masses  co  be  measured.  The 
Macfarlane  group,  for  example,  has  recently  detected  the  molecular  Ion  of  a  protein  with  m/e 
of  greater  than  12,000  (195). 
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Figure  60  -  Desorption  due  to  energetic  collision  cascades. 

Anocher  approach  has  been  che  use  of  an  aeon  beam  rather  than  an  Ion  beam  to  serve  as 
the  bombarding  species  (196,197).  Probably  for  commercial  reasons,  the  writers  have  termed 
the  technique  FAB  or  fast  atom  bombardment  rather  chan  SIMS.  The  physics  of  che  tw 
processes  are  completely  Identical  and  It  Is  not  clear  why  tw  names  are  needed  (See  Figure 
60).  On  che  ocher  hand,  the  FAB  source  offers  two  unique  practical  advantages  which  promise 
to  extend  che  SIMS  applications.  First,  the  use  of  an  atom  beam  allow  the  secondary  Ions  to 
be  extracted  from  che  sample  with  large  electrostatic  potential.  The  use  of  a  primary  Ion 
would  not  be  as  straightforward  since  It  would  be  deflected  by  such  a  field.  Secondly,  the 
proponents  of  che  FAB  source  suggest  that  for  insulating  samples,  charging  is  greatly 
minimized  since  the  Incident  beam  does  not  alter  che  sample  charge.  It  is  not  yet  clear  chat 
this  apparenc  advantage  cannot  be  offset  by  the  use  of  standard  charge  compensation  tools 
such  as  sample  biasing  and  secondary  electron  flooding.  It  is  of  further  interest,  however, 
chat  these  workers  have  discovered  certain  low  vapor  pressure  liquids  like  glycerol  which  act 
as  ideal  solvents  for  che  sample.  For  chls  case,  che  dilution  by  che  matrix  apparently 
minimizes  recombination  and  serves  to  bring  a  continually  fresh  surface  to  che  incident  beam, 
virtually  eliminating  any  sample  damage  problems.  A  typical  SIMS  spectrum  taken  using  a  FAB 
source  and  magnetic  analyzer  Is  shown  in  Figure  61 . 
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Plgure  61  -  Fast  atom  bombardment  spectrum  o£  undecapeptide,  methlonyl- 
lysl-brady-klnln,  obtained  from  l  ug  o£  material.  From 
reference  197. 
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